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ABSTRACT 

The formation of planets depends on the underlying protoplanetary disc structure, which in turn influences both the accretion and 
migration rates of embedded planets. The disc itself evolves on time scales of several Myr, during which both temperature and density 
profiles change as matter accretes onto the central star. Here we used a detailed model of an evolving disc to determine the growth 
of planets by pebble accretion and their migration through the disc. Cores that reach their pebble isolation mass accrete gas to finally 
form giant planets with extensive gas envelopes, while planets that do not reach pebble isolation mass are stranded as ice giants and 
ice planets containing only minor amounts of gas in their envelopes. Unlike earlier population synthesis models, our model works 
without any artificial reductions in migration speed and for protoplanetary discs with gas and dust column densities similar to those 
inferred from observations. We find that in our nominal disc model, the emergence of planetary embryos preferably tends to occur after 
approximately 2 Myr in order to not exclusively form gas giants, but also ice giants and smaller planets. The high pebble accretion 
rates ensure that critical core masses for gas accretion can be reached at all orbital distances. Gas giant planets nevertheless experience 
significant reduction in semi-major axes by migration. Considering instead planetesimal accretion for planetary growth, we show that 
formation time scales are too long to compete with the migration time scales and the dissipation time of the protoplanetary disc. All in 
all, we find that pebble accretion overcomes many of the challenges in the formation of ice and gas giants in evolving protoplanetary 
discs. 
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1. Introduction 

The formation of planets takes place in protoplanetary discs that 
surround newly born stars. This process can happen on differ¬ 
ent time scales. Gas giants have to form within the lifetime 
of the gaseous protoplanetary disc, since they must accrete a 
gaseous envelope after the formation of the planetary core. Be- 
eause the typical disc lifetimes are constrained to a few Myr 
iHartmann et al.lfT998t iHaisch et al.ll200lt lMamaiekll2009h . gi¬ 
ant planet formation has to happen within the same time span. 
Terrestrial planet formation, on the other hand, finishes on much 
longer time scales. In our own solar system the last giant impact 
is constr ained to have occurred ~ 100 Myr after t he solar system 
formed (iKleine et al.ll2009l:I.Tacobson et al.ll2014t) . 

The evolution of the protoplanetary disc is crucial for the for¬ 
mation of giant planets, because the whole growth process from 
dust to planetary cores happens in the gas phase of the disc. The 
structure of the disc can be approximated by a minimum mass 
solar nebular (MMSN), in which the gas surface density Eg, tem- 

S ierature T, and aspect ratio H/r are set to be simple power laws 
Weidenschillind fl 977t lHavashilll98lh . However, the structure 
of the inner regions of protoplanetary discs (r < 10 AU) is in 
reality much more complicated than a simple power law, and 
the disc ins tead features bumps and dips caused by transitions 
in op acity ( Bitsch et^l2014i iBitsch et~^ 120151: iBaillie et al.l 


in op acit y (IHitsch et al.M/Ul4t iHitsch et al.l l2Ul5t iHaillie et al.l 
12015li . In iBitsch et al. ( 2015li we ran an extensive suite of proto¬ 
planetary disc models and provided analytical fitting models of 
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the disc structure that greatly improve on the simplified MMSN 
model. 

We now review some critical processes for how planets are 
formed and how they are influenced by the structure of the pro¬ 
toplanetary disc. 

Planetesim al formation can be triggered through the stream¬ 
ing in stability (lYoudin & Goodman 200^ Llohansen & Youdinl 
l2()07h or in vortices (iRaettig et al.l 1201^ For this process to 
happen, dust particles in the protoplanetary disc f irst have to 
grow to pebbles by coagu l ation and condensation (IZsom et al.l 
I 2 OIOI: iBirnstiel et al.l 1201^ iRos & Johansenll2013l) . These peb- 
bles then move radially towards the sta rt owing to gas drag 
(IWeidenschillingll977HBrauer et al.l2008l) . During their motion, 
a swarm of pebbles can undergo a gravitational collapse and 
form a planetesimal, which are the first building blocks of plan¬ 
ets. The formation process of planetesimals via the streaming in¬ 
stability is a strong function of the pressure gradient in the disc, 
making discs with bumps in their pres sure profile very ap pealing 
compared to the simple MMSN disc (iBitsch et al.ll2015li . 

The growth of the core of a giant planet can happen via the 
accretion of planetesimals onto a planetary embryo, which is 
basically a large planetesimal. However, this process can eas¬ 
ily t ake longer than the lifetime of the protoplanetary disc it¬ 
self (iPollack et al.lll99^ lRafiko'^12004 [Levison et al.ll201(ili . if 
the amount of solids is not increased by a factor of 6 - 8 com¬ 
pared to solar value. These growth time scales can nevertheless 
be significantly reduced when the accreti on of small pebbles 
onto p lanetesimals is taken in t o account ( Johansen_&Lacer^ 
201TL_lOrmel_&Kl^ 201^ iLambrechts & JohansenI 2012 : 


Morbidelli & Nesvornvll2012 1. In this process the formation of 
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a core of ~ lOMearth at ~ 5 AU can occur within 1 Myr. Pebbles 
can either form by coagulation of dust in the disc or be the re¬ 
sult o f fragmenta tion of the planetesima l population (IChambersI 
I 2 OI 4 I 1 . Recently, Llohansen et al.l (1201 5h have also shown that 
even chondrules can be accreted effectively on planetesimal 
seeds. The process of pebble accretion depends signihcantly on 
the disc structure because the accret ion rate depends on the disc’s 
aspect ratio and pressure gradients (iBitsch et al.ll2015l) . 

Gas accretion starts after the core has reached its isola¬ 
tion mass for either planetesimals or pebbles, because during 
the accretion of solids the atmosphere is heated by the im¬ 
pacts, preventing a contraction and efficient gas accretion. This 
isolation mass depends o n the semi-major axis and column 
density of plan etesimals (iMizunol 119801: iKokubo & Idal 120021: 
[Raymond et al.ll2014l) and on the disc’s aspect ratio Hjr for peb- 
ble accretion (iLambrechts et al.ll201^ . After the bombardment 
of planetesimals or pebbles has stopped, the gase ous envelope 
can c ontract, and runaway gas accretion can start (iPollack et al.l 
Il996h . However, the isolation mass of planetesimals is very hard 
to determine and much clearer to determine for pebble accre¬ 
tion. For planetesimal accretion, t he isolation mass ha s to be 
estimated by N-body simulations (iKokubo & Idal 1200^ . while 
for pebble accretion the isolation mass is determined directly by 
the modifications of the pressure gradient in the disc, which is 
caused by the planet itself, and that can halt pebble acc retion 
(iPaardekooner & Mellemal I2OO6I: ILambrechts et al.l 120141) . The 
gas accretion is limited not only by the properties of the planet 
(e.g. mass of the core), but also by the disc itself, because the 
planet cannot accrete more gas than what is provided through 
the accretion rate of the gas onto the central star. In fact, the ac¬ 
cretion rate onto the plan et is at its maximum rough ly ~ 80% of 
the stellar accretion rate (iLubow & D’AngeI^l2006l) . 

Planet migration describes the gravi tational inte ractions of 
the planet with the surrounding gas disc (IWardlll99^ . In the lo¬ 
cally isothermal limit, the time scales for inward migration of 


embedded planets is shorter than the disc’s lifetime (t„ 


8 X 


10^ yr for an Earth size pla net at 5 AU), which poses a problem 
for the formation of planets (iTanaka et al.l2002ll . Considering the 
thermodynamics inside the dis c, recent studies have shown that 
planets can migrate outwards (Paardekooper & Mellema 


iklev & Cridall2008HBaruteau & Massetll2008UKlev et al.feoO' 

This outward migration depends on the gradient of entropy in 
the disc and is most likely to happen in regio ns of the disc wher e 
Hjr drops with increasing orbital distance (iBitsch et al.ll2013h . 
Low-mass planets are in this so-called type-I-migration phase, 
where the perturbation of the planet onto the disc is small. When 
the planets become more massive, for example, when it is caused 
by gas accretion, they start to open a gap inside the disc and mi¬ 
grate with the viscous accretion speed of the disc, which is much 
slower than the type-I rn igration and is called type-II migration 
(iLin & Papaloizoulll986l) . Migration thus affects planets of all 
masses, where its effects become significant when the planet is 
larger than one Earth mass. 

All these aspects and their interplay have to be considered 
when trying to explain the observed distribution of planets and 
exoplanets. Eirst attempts to explain the distribution of exoplan¬ 
ets have been done in so ca lled population synthesis studies tha t 
started about a decade ago (llda & Linll2004tlAlibert et al.ll2004ll . 
These studies generally assume that the core grows via the ac¬ 
cretion of planetesimals, after which gas accretion can set in. 
During this growth phase, the planets migrate through the disc. 
These models are able to explain the distribution of the observed 
exoplanets only by making some critical assumptions, some of 
which are questionable. These questionable assumptions regard¬ 


ing the migration speed of planets, the amounts of solids in the 
disc and the lifetime and evolution of protoplanetary disc it¬ 
self are discussed in more detail in section lb!9l where we show 
that no supposedly helpful assumptions have to be made in our 
model. 

The aim of this paper is to study the formation and evolu¬ 
tion of planets in evolving accretion discs around young stars, 
where planets first grow via pebble accretion and can then con¬ 
tract a gaseous envelope. We focus here on the formation of dif¬ 
ferent planetary types that can emerge in protoplanetary discs, 
on the parameters in initial semi-major axis, and on the initial 
tim e needed to form p la nets of a certai n planetary type. _ 

I^mbrechts et al.l (1201 4l) and ILambrechts & .lohansenl 
(l2()ll ~ propose that the dichotomy between ice and gas 
giants is a natural consequence of fast growth by pebbles 
and the existence of a pebbl e isolation mass. Additionally, 
ILambrechts & JohansenI (1201 4l) show that ice giants can over¬ 
come the type-I migration barrier by growing faster than they 
migrate. In this study we wish to investigate this concept in a 
more realistic disc than an MMSN, compared to their study 
improved planet migration and gas accre tion rates. 

We use the disc evolution model of IBitsch et al.l (1201 5h for 
solar-type stars. This disc model is a semi-analytical formula ht- 
ted to 2D radiation hydrodynamic simulations that feature stellar 
and viscous heating, as well as radiative cooling. It reproduces 
the dips and bumps in the disc profile caused by opacity tran¬ 
sitions and captures the disc evolution on a time scale of sev- 
eral Myr, which are l inked to observations of accretion discs 
(iHartmann et al.l ll 9981) . In this disc we implant planetesimals 
that accr ete pebbles, fo l lowin g the radial-drift-dominated ap- 
proach of iBrrnstiel et al.l (1201 2h ; ILambrechts & .lohansenl (1201^ 
for the formation of pebbles. These planets grow rapidly a nd 
can reach their pebble isolation mass (ILambrechts et al.l2014t) in 
several 100 kyr, which is when their gas accretion starts. The gas 
accret ion is modelled by using accretion rates ofjMachida_e^^ 
(201T) and envelope contraction rates following iPiso & Youdin 
(l2014h . During their growth, the planets migrate through the 
disc. We use the analytical torq ue formula for type-I migra¬ 
tion of IPaardekooner et al.l (1201 ih to mimic their motion in the 
disc. When planets become massive and start to open up a gap 
in the disc, they migra te with the viscous type-II migration 
(iLin & Papaloizoul 19861) . which is slower than type-I migration. 

Our work is structured as follows. In section|2]we explain the 
different methods used for pebble and gas accretion, for the disc 
evolution, and for planetary migration. We then present results of 
simulations where the planets grow via pebble and gas accretion 
while they migrate through the evolving disc (section [3]). The 
results obtained with pebble accretion are compared with simu¬ 
lations where the cores grow via planetesimal accretion in sec- 
tion|4] In section|5]we discuss the formation of the giant planets 
in our own solar system via pebble accretion. The many applica¬ 
tions of our planetary growth model are discussed in section |6] 
We hnally summarize in section Q 


2. Methods 

The methods used in this work are explained in much more detail 
in the literature cited in the following paragraphs. This section 
only intends to summarize the methods in a condensed way, so 
that it is easy to understand the principles on which our work 
is based. During the disc evolution in time (section IZTI) . planets 
grow first via pebble accretion (section |2^ very quickly. After 
they have reached their pebble isolation mass, gas can accrete 
onto the planet (section |23]l. During the whole growth process 
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planets migrate through the disc, which changes their semi major 
axes (section l24l) . 


2.1. Evolution and structure of the disc in time 


The lifetime of pro t oplanetary discs s pans 1 -10 My r 
(iHartmann et alJ 119981: iHaisch et alJ 120011: iMamaiefl l2009l) . 
During the lifetime of the disc, the accretion rate M changes 
in time following constraints from observati ons of slightly 
sub-s olar mass stars in the Taurus cluster (IHartmann et alJ 
Il998h . 


log 


M \ 
MJyrj 


= -8.00- 1.40 log 


t+ 10^yr\ 
lO^yr ) ■ 


( 1 ) 


The accretion rate M can then be related to the viscosity v and 
the gas surface density Zg via 


M — 3nvl,g , (2) 

where we assume a constant accretion rate for each or- 
bital distance. For the v iscosity, we take the a approach 
dShakura & Sunva^ll973h with a = 0.0054 constant through¬ 
out thewhol£disCj_w^ere V = 

iBitsch et'aP (12015l) calculated the structure of accretion discs 
around solar type stars with 2D simulations that includ viscous 
and stellar heating, as well as radiative cooling for several dif¬ 
ferent M rates, which correspond to dif ferent evolution tim es of 
the disc (Eq.[T]i. We note here that in the lBitsch et alJ(l2015l) disc 
model, the a value only represents the heating of the disc and 
is not representative of the viscous evolution of the disc. They 
then provided a semi-analyt ical fit to the disc str ucture evolution 
in time (see Appendix A in iBitsch et al.l (12015l) l. which we use 
for the evolution model of our disc. This model covers a radial 
extent from 1 to 50 AU. Inside of one, we extrapolate the fit of 
iBitsch et~aD (120151) with the given power laws. 

This extension of the disc structure fit is correct as long as 
the temperature in the disc is so low that silicates do not melt or 
evaporate. The melting or evaporation causes an additional tran¬ 
sition in the opacity profile, which changes the cooling prop¬ 
erties of the disc and therefore the structure of the disc. In the 
very early stages of the disc, the silicate evaporation line is at 
0.7 AU, but it moves inwards in time as the disc loses mass, so 
that silicates only evaporate in the very inner regions of the disc 
(r < 0.1 AU). In addition, we focus on planets that form in the 
outer disc (rp > 3 AU), which only reach the inner regions of 
the disc via migration when they have stopped accreting solids 
(Eg.lTbl) in the first place^ _ 

The disc structure of iBitsch et al.l (12015h features bumps and 
wiggles in the important disc quantities (Eg, T, and H), which 
are caused by transitions in the opacity k (e.g. at the ice line) that 
influences th e cooling rates of the disc as D oc 1/a: (IBitsch et alJ 
1201 31 l2014l) . A change in the cooling rate of the disc directly 
changes the discs temperature T and thus the scale height of the 
disc {T oc (///r)^], which in turn changes the local viscosity of 
the disc. This change in the local viscosity has to be compensated 
for by a change in the surface density Eg to have the same M at 
all orbital distances, thus creating a chang e in the local radia l 
gradient in surface density and pressure P (IBitsch et al.ll2014l) . 
Therefore a steeper gradient in temperature will result in a shal¬ 
lower gradient in surface density at the same orbital location. 

This has important consequences for the accretion of pebbles 
(section lZ2l i. which depends on the pressure gradient parameter 
T] (eq. |6]l and for the migration of planets, which depends on 


the gradients of temperature, surface density, and entropy (sec¬ 
tion | 23 ]). 

The disc structure in itself depends on the dust grains inside 
the discs, because those grains are responsible for the absorp¬ 
tion and re-emission of photons that distribute the heat inside 
the disc. The main contribution to the dust opacities originates in 
micrometre-sized dust grains. Larger dust grains only contribute 
minimally to the opacity, so that we do not take their contribution 
into account. We assume here a metallicity of micrometre-sized 
dust grains of 0.5% or 0.1% of the gas density at all time. We also 
make the assumption that this small dust is coupled perfectly to 
the gas and does not ev olve its size distrib ution in time. Here we 
use the opacity table of iBell & LinI(Il994l) . 

The decay of the disc accretion rate from M -lx lO^^Mp/yr 
down to M - \ X 10 ®Mo/yr takes 5 Myr in IHartmann et al.l 
(Il998h . Using the time evolution of M via eq. [T] the disc spends 
2 Myr decaying from M - 2x 10^®Mo/yr to M - lx 10“®Mo/yr. 
However, for these low accretion rates, photoevaporation be¬ 
comes very efficient, and the disc c an dissipate in much shorter 
time scales (lAlexander et al.ll2014l) . Eor this reason, our nomi¬ 
nal disc lifetime is set to 3 Myr, which is when we assume that 
photoevaporation clears the disc immediately, but we follow the 
decay rate of M given by eq.[T]down to M -2x 10“®Mo/yr. 

The disc structure significantly changes as the disc evolves 
in time and as M decreases. As the disc reduces in M and Eg, the 
disc becomes colder, because viscous heating decreases, which 
implies that the opacity transition at the ice line moves inwards. 
This means that the bumps and wiggles in the disc structure 
[T, Eg and H) move inwards as well. The star also evolves and 
changes its luminosity, changing the amount of stellar heating 
received by the disc and thus changing the temp erature, which 
is all taken into account in the lBitsch et al.l (1201 5l) model. These 
changes to the disc structure influence the formation and migra¬ 
tion (see Eig. [T]i of growing protoplanets significantly. 


2.2. Growth via pebbles 


T he growth of planetary embryos v ia pe bble accretion is outlined 


in 


___ Lambrechts & JohansenI (12012l) and iLambrechts & JohansenI 
(l2014h . The pebbles form from grains initially e mbedded in the 
proto planetary disc (~ qm size) by collisions (iBirnstiel et al.l 
12012 l) or through sublimation a nd condensation cycles around 
ice lines (iRos & Johansenll2013l) . Swarms of these pebbles drift 
inwards towards the star, but can collapse under their own grav¬ 
ity and form planetesimals of 100 to 1000 km in size in a 
process called streaming i nstability (lYoudin & Goodm^l2005l : 
Llohansen & Youdinll2007l) . E urther discussion on this process 
can be found in the review of iJohansen et al.l (l2014l) . along with 
a list of other models of planetesimal formation by particle con¬ 
centration and gravitational collapse. 

We now consider cores that predominantly grow by accre¬ 
tion of particles with approximately mm-cm sizes. This particle 
size can be expressed through the gas drag time scale tf and the 
Keplerian frequency Qk in terms of the Stokes number 


Tf = PlKtf = 


p.R 

Pg^g 


(3) 


where p, is the solid density, R the particle radius, pg the gas 
density, Qk the Keplerian frequency, and Hg the local gas scale 
height. Small particles (Tf <s; 1) are strongly coupled and move 
with the gas, while larger particles (Tf » 1) are only weakly 
affected by gas drag. 

The scale height of pebbles //peb is related to the scale height 
of the gas Hg through the viscosity and the Stokes number 
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dYoudin & Lithwicfl2007l) by 

Hpeh = Hg ^Ja/Tf, (4) 


where a is the viscosity parameter. In our simulations we place 
seed masses that have reached the pebble transition mass 


Ml = 


(r]VK)^ 

V 3 GQk ’ 


where G is the gravitational constant, uk = ^Kf, and 

1 ///\^ dlnF 
^ ~ ”2 17/ d\nr ■ 


(5) 


(6) 


Here, 5 In P/5 In r is the radial pressure gradient in the disc. 
These masses are typically in the range of 5 x IO^^Me to IO^^Me 
( see Fig. |2|i. The pebble transition mass defines the planetary 
mass at which pebble accretion occurs within the Hill radius, 
while for lower masses pebbles are a ccreted within the Bondi 
radius (iLambrechts & .lohansenllTlT^ . 

These masses are a bit higher than planetesimals formed 
by the streaming in stability, which have roughly 10 “"^ Me 
( iJohansen et al.llTlT^ . Even if a planetesimal of 10 Me forms 
at f = 0, the planetesimal has several Myr to grow to the peb¬ 
ble transition mass Mi. This growth phase can occur through 
the accretion of p lanetesimals or pebbles in the inefficient Bondi 
accre tion regime (ILambrechts & JohansenllTlT^ Llohansen et al.l 

IMl. 

Planets whose Hill radius is roughly larger than the scale 
height of the pebbles (eq. nil accrete in a 2D fashion, and the 
accretion is given by 

^c,2D = '■H^H^peb , (7) 

where th = r[Mc/(3M*)]*^^ is the Hill radius, uh = the 
Hill speed, and Epeb the pebble surface density. If the Stokes 
number of the particles Tf is larger than 0.1, the accretion rate 
is limited to 


4^fc,2D - 2rHUH2ipeb , 


(8) 


because the pla netary seed cannot accrete part icles from outside 
its Hill radius (ILambrechts & Johansenll2012h . However, when 
the planets are small and their Hill radius is smaller than the scale 
height of the pebbles, the pebble accretion rate is reduced and 
plan ets accrete in a 3D way , which is related to the 2D accretion 
rate (iMorbidelli et al.ll2015h by 


^c,3D 


- ^c,2D 


;r(Tf/0.1)‘^VH ' 
, 2 VTrFfpeb . 


(9) 


The transition from 3D to 2D pebble accretion is then reached 
(IMorbidelli et al.ll2015[) when 


7r(Tf/0.1)‘^^rH 

2y/2^ 


^ 7/peb . 


( 10 ) 


This transition depends on particle size (rf) and on the scale 
height of the disc. This means that in the outer parts of the disc, 
where //peb is larger and rf is smaller, a higher planetary mass is 
needed to reach the faster 2D pebble accretion branch. We use 
the Stokes number of the dominant particle size 


V3 tp Lpeb 

8 ?7 Eg ■ 


( 11 ) 
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This size is obtained from an equilibrium between growth 
and drift to fit constraints from advanced coagulation mod- 
els and observation s of pebbles in protoplanetary discs 
( Birnstiel et al. 7) T^. The par ameter ep is 0.5 and eo is 0.05 


( Lambrechts & Johansenl2()14l) . In our disc model, this results in 
Stokes numbers between 0.05 and 0.5. The pebble surface den¬ 
sity depends on the gas surface density Lg and the semi major 
axis Tp of the planet through 


2MpebLg 

y V37r€prpVK 
where the pebble flux is 
dr„ 

Mpeb = 27rrg-^(ZLg) . 


( 12 ) 


(13) 


Here, Z denotes the fraction of solids (metallicity) in the disc that 
can be transformed into pebbles at the pebble production line rg 
at time t 

r, - (GM,yi\euZflh'^l^ , (14) 


$ " ^ , (15) 


where M* is the stellar mass, which we set to IMq. After 3 Myr 
of disc evolution, the pebble production line is at 250 AU, in¬ 
dicating that the disc has to be at least 250 AU wide to sus¬ 
tain a pebble flux for 3 Myr. After 5 Myr of disc evolution, the 
pebble produc t ion lin e is located at 360 AU. Observations by 
lAndrews et al.l (1201 Ol) And typical protoplanetary disc radii to be 
150 - 200 AU for discs that are a few Myr old, which is only 
slightly smaller than our estimated pebble production line at 3 
Myr. 

When pebbles form in the outer disc drift across the wa¬ 
ter ice line, they will melt, because they mainly consist of ice, 
and release the trapped silicate particles. Then the particle size 
shrinks significantly, which will slow down the pebble accretion 
rate onto the planet (eq. |7]i. However, in our model, the water 
ice line is located outside 3 AU only in the very early stages 
of th e disc evolution, a nd it moves inwards very quickly with 
time (iBitsch et al.l2015h . Since planets grow rapidly by accretion 
of pebbles locally and experience most migration after reaching 
pebble isolation mass, the planets only accrete icy pebbles and 
reach their pebble isolation mass before they migrate across the 
water ice line. 

The initial planetary seeds efficiently accrete pebbles (eq. |2l) 
and grow very fast. During this growth process, the planetary 
seeds can also attract a gaseous envelope (section l2jT l. The struc¬ 
ture of the gaseous envelope is supported by the accretion lumi¬ 
nosity deposited by the accreted pebbles into the atmosphere of 
the planet. The envelope will collapse when the mass of the core 
is similar to the mass of the envelope itself. However, when the 
envelope collapses, this critical core mass is a function of the 
accretion rate onto the planet itself, where the critical core mass 
becomes higher with increasing accretion rates. Without any in¬ 
terruption in the pebble accretion rate, the critical core mass can 
be up to ~ IOOMe, which is up to an order of magnitude higher 
than the amount of solids in the giant planet’s cores of the solar 
system (ILambrechts et al.ll2014t) . 
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However, when the planet reaches a certain mass, it changes 
the gas pressure gradient in the disc locally, which modifies the 
rotation velocity of the gas that then halts the drift of pebbles that 
could be accreted onto the core so the accretion of pebble s stops 
(iPaardekooper & Mellemal [20061: iLarn brechts et alj 201 ^. Th is 
is the so-called pebble isolation mass ( Lambrechts et al.ll2014h 


Mi,o^20{^) MEarth. 


(16) 


When pebble isolation mass is reached, the planet can contract 
its envelope and start gas accretion. The pebble isolation mass 
is therefore a natural division between gas and ice giants, where 
ice giants did not reach pebble isolation mass early in the disc 
lifetime. This can then explain th e difference between ga s and 
ice giants in our own solar system (ILarnbrechts et al.ll2014ll . 


2.3. Gas accretion 


Planets that reach pebble isolation mass (eq. [Thl l can start to ac¬ 
crete gas, because the pebble flux onto the planet that heats the 
envelope and hinders its contraction has stopped. However, dur¬ 
ing the formation of the planetary cores via pebble accretion, 
small amounts of highly polluted gas can be bound to the planet 
inside the planets Hill sphere. We therefore assume that 10% 
of the planets mass is in gas, prior to the point when the planet 
reaches pebble isolation mass. This means that 90% of the nomi¬ 
nal pebble accretion rate is counted as solids, and 10% is counted 
as gas. During this stage the planet grows at the same speed as 
if 100% of the accretion were in pebbles. This approach also 
means that the same pebble isolation mass is reached, because 
the total mass of the plane t is the same. 

iPiso & YoudinI (1201 4l) estimate the gas contraction time of a 
gaseous envelope around a planet. After the planet has reached 
pebble isolation mass, the envelope of the planet contracts on 
a long time-scale while it accretes some gas. This contraction 
phase takes place as long as Menv < and the corresponding 
gas accretion rate can be extracted from iPiso & YoudinI (l2014ll 
and is given by 


(O.IMe) (SIK; Myr’ 


(17) 


where / is a fudge factor to change the accretion rate in order 
to match nu merical and analytical results, which is normally set 
to / = 0.2 (iPiso & YoudinluOf^ . The opacity in the planets 
envelope /Cenv is generally very hard to determine because it de¬ 
pends on the grain sizes and distribution inside the planetary at¬ 
mosphere. Here we use /Cenv = 0.05cm ^/g, which is very sim¬ 
ilar to the values used in the study by iMovshovitz & PodolakI 
(l2008l) . In Appendix iBl we test the influence of different values 
of /Cenv for gas accretion. For the density of the core, we as¬ 
sume Pc = 5.5g/cm^. This contraction phase ends as soon as 
iHcore = iHenv and rapid gas accretion starts. 

For rapid gas accretion (Mcore < we follow 

iMachida et al.l (1201 Oh directly. They calculated the gas accretion 
rate using 3D hydrodynamical simulations with nested grids. 
They find two different gas accretion branches, which are given 
as 

Mgns.low = O.SSflKSg//" (18) 

and 

Mgas.high = 0.14QK£g//^ (19) 


where the effective accretion rate is given by the minimum of 
these two accretion rates. The low branch is for low mass plan¬ 
ets (with (RhIH < 0.3), while the high branch is for high mass 
planets {{Rulh > 0.3), and the effective accretion rate is given by 
the minimum value of both rates. Additionally, we limit the max¬ 
imum accretion rate to 80% of the disc’s accretion rate onto the 
star, be cause gas can flow through the gap, even for high mass 
planets (iLubow & D’Ang^l200^ . 


2.4. Planet migration 

The growing protoplanets inside the disc interact with the sur¬ 
rounding gas and migrate through it. The process of migration 
is substantially different between low mass planets that are still 
fully embedded in the disc (type-I migration) and high mass 
planets that open up a gap inside the disc (type-II migration). 
The migration rates of low mass planets can be obtained by 
2D and 30 hydrodynamical simulations (iKlev & Cridal 120081: 
iKlev et aklllKM iBitsch & Kle^l201 ll: iLega et al.ll2()T4ir How¬ 
ever, these simulations are very computationally intensive, so we 
use a prescri bed formula to compute the torque acting on embed- 
ded planets dPaardekooperet al.ll2()Tlh . The torque formula of 
IPaardekooper et a" (201 ih includes the effects of torque satura¬ 
tion a nd has been tested ag ainst 3D simulations in fully radiative 
discs (lBitsch_&Klev|^0^, which find good agreement. Recent 
studies of iLega et ^ (120151) tested the torque formula against nu¬ 
merical simulations in accreting discs including stellar and vis- 
cous heating and r adiative cooling, as used in the disc model of 
iBitsch et al.l (l2015l) . and found very good agreement with respect 
to the zero-torque location in the disc, where planets would stop 
their inward migrat ion. For low mass planets (Mp < SMearth), 
previous studies of iLega et al.l (l2014t) have shown a slight dis - 
crepancy with the torque formula of IPaardekooper et alJ (1201 iT) . 
however these differences were found to be very small, even con¬ 
sidering that very small planets migrate very slowly in the first 
place. Here we also assume that planets move only on circu¬ 
lar orbits around the stars, because eccen tricity and inclination 
is damped quite qui ckly by the gas disc (iBitsch & Kiev! 120101 : 
IBitsch & Kle\ll20]Th . 

We therefore use the torque formula of IPaardekooper et al.l 
(1201 ll) . where the total torque Ftot acting on a planet is given as a 
composition of the Lindblad torque Fl and the corotation torque 

Tc, 


Ttot = Tl + Fc . (20) 

The Lindblad and corotation torques depend on the local radial 
gradients of surface density Lg oc temperature T oc and 
entropy S oc with ^ - p - {y - L0)s, where y = 1.4 is the 
adiabatic index. 

Very roughly said, for 2g gradients that are not too negative, 
a radially strong negative gradient in entropy, caused by a large 
negative gradient in temperature (large /?), will lead to outward 
migration, while a shallow gradient in entropy will not lead to 
outward migration and planets migrate inwards. Therefore plan¬ 
ets can migrate outwards in certain regions of the disc, where 
strong negative gradients in temperature can be found. Generally 
these regions of outward migration exist close to transitions in 
opaci ty, where H/r drops (IBitsch et al.ll2013[ 120141: IBitsch et"^ 
l2015l) . However, as the disc evolves in time, these regions of out¬ 
ward migration also evolve in time, so that at the very late stages 
of the disc evolution (M < 4 x 10“®Mo/yr) only very small re¬ 
gions of outward migration exist that can only hold planets of up 
to ^ lOMEarth- 
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1I2OI ih for our standard disc structure with Z^ust = 0.5%. If the planet 


is located inside the solid lines, it will migrate outwards. If the planet 
is outside the solid lines, it will migrate inwards. The small ticks on the 
top of the plot correspond to the location of the ice line at the times 
indicated by the different colours. 


This is illustrated in Fig. [T] where we display the evolution 
of the regions of outward migration in time in our simulations. 
As the disc evolves in time, the strong negative gradient in en¬ 
tropy, which is caused by the negative gradient in temperature 
just outside the ice line (r > rice) that can trigger outward mi¬ 
gration, moves towards the star. The region of outward migra¬ 
tion therefore moves towards the star, which finally makes the 
inner region of outward migration (caused by the silicate line) 
disappear after 0.75 Myr. In the late stages outward migration 
is only possible in the inner parts of the disc for very low mass 
planets (Mp < lOMp), where the regions of outward migration 
stay roughly constant, because the disc’s temperature gradients 
do no t evolve significant ly at these stages of disc evolution any 
more (iBitsch et al.ll2015T) . 

Planets that have reached their pebble isolation mass leQ.fTbll 
start to accrete gas (section l23T l and grow even further until they 
finally open a gap inside the disc. A gap can be opened (with 
^icap < 0.12g), when 


^ 3 // 50 , 

V =-H- — < 1 , 

4rH cfR 


( 21 ) 


where rn is the Hill radius, q — MpjMj,, and 7? the Reynolds 
number given by 7? = rpf2p/v dCrida et al.ll2006l) . If the planet 
becomes massive enough to fulfil this criterion, it opens up a 
gap in the disc, and it migrates in the type-II regime. The gap¬ 
opening process splits the disc in two parts, which both repel the 
planet towards the centre of the gap, meaning that the migration 
time scale of the planet is the accretion time scale of the disc, 
Tvisc = r-pjv. However, if the planet is much more massive than 
the gas outside the gap, it will slow down the viscous accretion. 
This happens if Mp > 47rSgrp, which leads to the migration time 
scale of 


/ 

Tii = Tv X max 1, 
V 


Mp ' 
47rSgr2, 


( 22 ) 


resulting in slower inward migration for massive planets 
(iBaruteau et al.ll2()T4l) . 


Before the planet is massive enough to open up a deep gap 
inside the disc, there is still material inside the planets corota¬ 
tion region, which reduces the total negative torque acting on the 
planet and can in principle be so strong t hat the planet can mi¬ 
grate outwards ( Crida & Morbidenill2007h . The depth of the gap 
is given in ICrida & Morbidellil ( 200% as 


/(n = 


P-0.541 



if !P< 2.4646 
, otherwise. 


(23) 


The factor f{'P) reduces the migration rate of the planet, when 
a partial gap is opened in the disc, because the migration rate 
depends directly on the gas surface density Sg. We multiply our 
migration rate directly by fiV) to reduce the migration rate by 
the partial opening of the gap. This reduction of the migration 
rate is very crucial, because planets that open a partial gap in 
the disc generally have several 10s of Earth masses and migrate 
very fast through the disc, because they are still in the type-I 
migration regime. 

In additional, we use a linear smoothing function for the tran¬ 
sition between planets that open partial gaps inside the disc (that 
migrate with the reduced type-I speed by the factor /(!P)) and 
planets that migrate with type-II, because even the reduced type- 
I migration rate (if the gap is fully opened with V < 1) is different 
from the nominal typ e-II rate. _ _ 

A new study by iBenltez-Llambav et al.l (1201 5h shows that 
low mass planets (Mp < 5 Mp) that accrete very fast can ac¬ 
tually migrate outwards instead of the inward type-I-migration. 
However, we find that this effect is not that important for plan¬ 
ets that grow via pebble accretion, because the growth is so fast 
that 5 Me is reached in a very short time, which also limits the 
time the planet actually migrates until it reaches 5 Mp. This is 
discussed in Appendix iDl 

The corotation torque arises from material that executes 
horseshoe U-turns relative to the planet, where most of this mate¬ 
rial is trapped in the planet’s horseshoe region. But, if the planet 
migrates with respect to the disc, material outside the horseshoe 
region will execute a unique horseshoe U-turn relative to the 
planet, which can alter its migration speed. This becomes impor¬ 
tant, in particular, when the planet starts to carve a gap around 
its orbit. This can lead to runaway type-III migration, if the coor¬ 
bital mass deficit 6M is greater than the mass of the planet Mp, 
which can significa ntly change the semi-major a xis of the planet 
in just a few orbits (iMasset & Papaloizoull2003h . The co-orbital 
mass deficit is defined as the mass that the planet pushed away 
from its orbit compared to the unperturbed disc structure as it 
starts to open up a gap. Unfortunately, there are no prescrip¬ 
tions to model this migration analytically, but we nevertheless 
test whether a growing planet might be subject to runaway type- 
III migration during our simulations. 

Migrating planets also experience dynamical torques that 
are proportional to the mig ration rate and depen d on the back¬ 
ground vortensitjQ gradient (lPaardekooperl20 1^ . These dynam¬ 
ical torques can have either positive or negative feedback on the 
migration, depending on whether the planet migrates with or 
against the direction of the static corotation torque. The effects 
of these dynamical corotation torques can be profound because 
they can slow down inward migration significantly, and outward 
migration can proceed beyond the zero-torque lines in discs that 
are massive and have a low viscosity. An approximate estimate 
of whether dynamical corotation torques play a role depends on 
the disc’s viscosity and mass, as well as on the planet’s mass 

* Vortensity is defined as the ratio of vorticity and surface density. 
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iPaardekooDej (l20 1 4l) . In particular, when 

vp.o 16/3 \ rtot/rog'jXj 

oQp,o ^ ^ I 2 “ 7 {Hjrf ’ 

migration due to the dynamical corotation torques will be impor¬ 
tant. Here the subscripts P and 0 indicate the initial location of 
the planet, s is the radial gradient of the surface density profile. 
To the normalisation of the torque, - 7rrpQZp,o/M* the disc 

mass, and - ^[qJh, where q - Unfortunately this 

effect is not quantified further, so that we just indicate when the 
dynamical corotation torques might become important, but our 
simulations do not evolve with them. This effect and the effect 
of type-III migration are both discussed in Appendix]^ 

As soon as the planet reaches an inner edge of 0.1 AU, we not 
only stop migration, but also the total simulation. That close to 
the central star, stellar tides can become important and influence 
the evolution of the planet, which we do not take into account. 


3. Formation of planets 

In this section we explore the growth of planetesimals that are 
inserted at a given initial time fp a given initial distance rp 
into the disc. The initial time to is important, because we keep 
the lifetime of the disc at 3 Myr, meaning that planets that are 
inserted into the disc, for example at fo = 2 Myr, will only expe¬ 
rience 1 Myr of evolution. Additionally, the structure of the pro - 
toplanetary disc is different at different ages dBitsch et alJl2015l) . 
The initial distance ro determines where the planetary seed is 
placed. This strongly influences the initial growth of the planet, 
because the surface density Sg is lower at greater orbital dis¬ 
tances, which means that the pebble surface density Speb is also 
lower, indicating a longer growth time of the core. After the plan¬ 
ets have reached their pebble isolation mass, the contraction of 
the gaseous envelope begins until Me < Menv At this point, run¬ 
away gas accretion starts. 

We focus here on discs that have a total lifetime of 3 Myr. A 
longer disc lifetime (5 Myr) does not affect our results qualita¬ 
tively, but simply pushes the preferred planet formation time out 
to ~ 3 Myr (see Appendix lAli. Additionally, in this section we 
assume that the opacity in the envelope is fixed to /Cenv = 0.05 
cm^/g. Different opacities of the envelope are explored in Ap¬ 
pendix |B] 

In Table [T] we define the different planetary categories used 
in this work. Our definition of different planetary categories is 
simply a function of planetary mass. Only the subcategories are 
a function of the final orbital d istance. The defin i tion o f ice gi¬ 
ants is slightly different than in lLambrechts et al] (l2014l) . where 
ice giants are required to not reach pebble isolation mass. This 
change in the definition is related to the slow contraction of the 
gaseous envelope (eq. [m, which allows for Me > Menv for a 
long time after reaching pebble isolation mass. 

3.1. Single evolution track 

In this section we follow the evolution tracks of planets in the 
disc to clarify the different outcomes of planetary evolution. We 
use a metallicity of Z = 1.0% in pebbles that can be accreted 
onto the initial seed masses. When we include the 0.5% of met¬ 
als in dust grains from the disc structure, our total metallicity is 
therefore roughly the solar value. We chose different initial start¬ 
ing locations ro of the planets, and all planets start at the initial 
time to = 2 Myr, which means that the planets will evolve for 1 



Fig. 2. Growth tracks of planets that accrete pebbles in an evolving pro¬ 
toplanetary disc. The big black circular symbols indicate the final mass 
and position of the planets at ff = 3 Myr, meaning that the planets ex¬ 
perience 1 Myr of evolution. The small black dots indicate 2.2, 2.4, 2.6, 
and 2.8 Myr. In the example of the purple line (ro = 10 AU), the planet 
reaches 0.1 AU before the final disc lifetime, and we stop the simula¬ 
tion. The different coloured lines indicate different planetary evolutions, 
which correspond to different types of planets, where we have an ice gi¬ 
ant, a hot gas giant, a cold gas giant, one more ice giant, and two ice 
planets starting from inside out (see text). The grey line indicates the 
surface density of pebbles in the disc at 2 Myr and ends at the ice line 
slightly interior to 1 AU, where the pebble surface density changes. . 


Myr, because our total disc lifetime is 3 Myr. In Fig.|2]the evolu¬ 
tion tracks of the planets and the pebble surface density Speb at 2 
Myr are shown. The surface density of pebbles Speb is calculated 
through eq. [12] which depends on Therefore the bu mps in 
2g caused by the transitions in opacity (iBitsch et ^12015h trans¬ 
late into bumps in the pebble surface density. As the gas surface 
density decreases and evolves in time, so does the pebble surface 
density. The starting mass of the planets is different at different 
locations, because it is set by the pebble transition mass, which 
depends on the disc’s aspect ratio (eq.|5l). 

Planets that manage to reach their pebble isolation mass con¬ 
tract their envelope until Me < Menv and the planets can start run¬ 
away gas accretion. During the whole growth process the planet 
migrates through the disc. The planets have been inserted into 
the disc at an disc age of fo = 2 Myr and are evolved until If = 3 
Myr. Planets that cross the ice line have already reached pebble 
isolation mass, so they do not accrete pebbles any more, making 
it unnecessary to model a transition in pebble size and pebble 
surface density at the ice line. 

The planetary seed implanted at 50 AU (green line in Fig.|2) 
accretes pebbles very slowly, because the pebble column density 
is very low at those orbital distances. At 3 Myr, the planet has 
just reached « 0.2 Me and is very far away from reaching pebble 
isolation mass. Because of its low mass, it also migrates only a 
very short distance. This is a typical example of an ice planet in 
our definition. 

The planet starting a bit farther inside at ro - 40 AU (light 
blue line in Fig. |2]i grows also very slowly, because the surface 
density in pebbles and the Stokes number (rf < 0.1) is low in 
the outer disc. It therefore never reaches pebble isolation mass 
and stays low, with a mass of only a few tenths of an Earth mass. 
Because the planet stays quite small, it only migrates a few AU 
inwards. This is an example of a slightly more massive ice planet 
in our definition, because it forms in the cold outer parts of the 
disc and is less massive than 2 Me. The planets starting at ro = 
40 and ro - 50 AU always accrete in the 3D scheme, because 
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Planet category 

Planetary mass 

Orbital distance 

Ice planet 

O.IMe < Mp < 2 Me 

- 

Ice giant 

Mp > 2 Me and M^ > Menv 

- 

Cold gas giant 

Me < Menv 

rf > 1.0 AU 

Warm gas giant 

Me < Menv 

0.1 AU<rf < 1.0 AU 

Hot gas giant 

Me < Menv 

rf < 0.1 AU 


Table 1. Definition of planets used in this work. We distinguish the main classes of planets (gas giants, ice giants, and ice planets) only by mass 
and mass ratio between core and envelope. The subcategories for gas giants are defined through their orbital distances, which we do not apply for 
ice giants and ice planets. 


(i) the particle scale height is large in the outer disc and (ii) Tf is 
very small, which explains their low growth rates. 

The planetary seed starting at 25 AU (black line in Fig. [2]i 
accretes pebbles at a faster rate (higher Epeb, larger Tf), but it 
misses reaching pebble isolation mass after 1 Myr of evolution. 
Because it is more massive, it migrates farther in the disc and 
ends up at ~ 22 AU. Its mass ends up at about 5 Me, where 
the mass of the core is much more massive than the mass of the 
envelope, making it an ice giant according to our definition. 

Starting a planetary seed at 15 AU (dark blue line in Fig.|2]i 
reveals a new growth path. After the planet has reached pebble 
isolation mass, its envelope contracts and the planet starts a run¬ 
away gas accretion process, making it a gas giant planet. Dur¬ 
ing its growth the planet migrates inwards from 15 AU down to 
~ 3 AU. The main migration happens when the planet is under¬ 
going fast inward type-I migration, before it is massive enough 
to open a gap in the disc (at ~ 120 Me). This fast inward migra¬ 
tion before gap opening indicates that the formation of gas giants 
requires a formation much farther out in the disc than their final 
orbital position would indicate; in situ formation of gas giants is 
impossible when including full planetary migration rates. 

The planetary seed starting at 10 AU (purple line in Fig. |2]i 
follows a similar evolution to the one starting at 15 AU, but 
with two notable exceptions. After the planet has reached pebble 
isolation mass, it migrates inwards very rapidly to 3 AU, with¬ 
out growing too much. This fast inward migration is caused by 
the disc structure, where the region between 3 and 8 AU has 
a very shallow (an d even inverted) radial temperature gradient 
(iBitsch et al.ll2015h . which causes a strong negative total torque 
acting on the planet driving fast inward migration. The planet 
is then in a region of very slow inward migration (just a bit too 
massive to be caught in the region of outward migration, Fig.[Tli, 
where it migrates very slowly and starts to rapidly accrete gas. 
During this accretion process, the planet migrates further into the 
inner regions of the disc. In fact, the planet migrates so fast that 
it reaches the inner edge of the disc at 0.1 AU, before the end of 
the lifetime of the disc is reached. (That is also why there is no 
big black circle in Fig.l^for this planet.) As soon as the planet 
reaches 0.1 AU, we stop the simulation. The planet has become 
a hot gas giant. 

If the planet starts in the very inner regions of the disc at 
5 AU (red line in Fig. [2]i, its isolation mass is very low, because 
H/r is very small in the inner parts of the disc at an evolution 
stage of 2 Myr. It therefore accumulates only 2.5 Earth masses 
of solids. A low core mass then leads to a very long contraction 
time of the envelope (eq. fTTl) . resulting in a total planet mass of 
only a few Earth masses when the disc reaches an age of 3 Myr. 
During the evolution, the planet is small enough to be caught at 
the zero-migration distance for most of its evolution, meaning 
that the planet follows the zero-migration distance as the disc 
accretes onto the star (see Fig. [1]). This planet is also classified 
as an ice giant because it formed at rp > rice (Fig- El- 


We would like to point out here that definitions of planets 
in this work is only related to the planetary mass and the mass 
ratio between the planetary core and envelope. The final orbital 
distance only plays a role in the subcategories of planets, so that 
gas giants very close to the central star (rf <0.1 AU) are called 
hot gas giants (TableE- 

3.2. Variation in the initial orbital position 

We can now expand Fig. |2]over the whole radial domain of the 
disc, but keep the initial time when we insert the planet fixed at 
k) - 2 Myr. At each orbital distance 3 AU < ro < 50 AU, a 
single planet is put into the disc and evolved independently. The 
nominal lifetime of the disc is 3 Myr, so the planets evolve for 1 
Myr in the disc. 

In Fig. [3] we present the final planetary core mass Me and 
envelope mass Menv, as well as total planetary mass Mtot with re¬ 
spect to the initial orbital position tq and the final orbital distance 
Tf of planets that evolved in the disc for 1 Myr after insertion at 
to -2 Myr. Planets inside the grey area are within rf <0.1 AU to 
the host star, and their evolution is stopped before 1 Myr of evo¬ 
lution is reached. The planetary masses and final positions are 
displayed as a function of their initial orbital distance tq, which 
illustrates what influence the initial orbital distance of the planet 
has on the final properties of the planet. For example, when look¬ 
ing at To = 15 AU, the total planetary mass Mtot and final orbital 
distance rf shown correspond to the location of the big black 
circle of the tq = 15 AU planet in Fig.|2] 

Planets forming outside of tq > 22 AU do not reach pebble 
isolation mass, because the pebble density in the outer disc is 
very low and growth time too long. But, these planets can still 
have a few Earth masses, making them ice giants by our defi¬ 
nition, because Me > Menv and Mp > 2 Me. However, planets 
starting with tq > 28 AU grow only very little, because of a 
small Tf and a larger scale height, which only allows planets to 
grow with the 3D accretion mechanism. These planets then have 
Mp < 2 Me and are ice planets instead. 

Just inside (20 AU < tq < 22 AU), the planets have reached 
pebble isolation mass, but did not reach runaway gas accretion, 
because Me > Menv, making these planets ice giants. Planets that 
form farther inside with 20 AU > tq > 13 AU become gas giants 
that stay in the outer disc with Tf > 1 AU. This indicates that 
there is a very broad range of radial extent that allows for the 
formation of gas giants, because their cores still form quickly 
enough via pebble accretion that enough time is left to accrete 
a gaseous envelope before the disc dissipates. With increasing 
To, the mass of the core increases, which is caused by the flar¬ 
ing st ructure of the disc at that evolutionary stage (iBitsch et al.l 
l2015h . which increases the pebble isolation mass (eq.fTbll. 

Planets that form within tq < 13 AU and outside of 6 AU 
all end up within Tf < 1 AU and are gas giants. The reason for 
their strong inward migration lies in the efficient growth of the 


Article number, page 8 of[23] 












Bitsch et al.: The growth of planets hy pehble accretion in evolving protoplanetary discs 



Fig. 3. Final planetary core mass iWcorc and corresponding envelope 
mass Mem, as well as total mass M^t = and final orbital 

distance rf of planets in the evolving disc at a disc dissipation time of 3 
Myr. The light blue line represents the initial orbital distance ro. Plan¬ 
ets inside the grey area are subject to very strong migration, so they 
end up in the inner parts of the disc with rf < 0.1 AU and become hot 
gas giants. The divide between ice and gas giants can be seen easily at 
~ 22 AU, where Me < Menv All data lines for one ro location in this 
plot give the final result of an individual growth track as in Fig.|^ 


planetary core via pebbles. Because the core grows very quickly, 
the planets have a longer time to migrate faster compared than 
cores that grow more slowly, because the type-I migration speed 
is proportional to the planetary mass. The planets then accrete 
gas to become gas giants and open up a gap in the disc that slows 
their migration (type-II migration). However, planets that form 
in a region of 7 AU< ro < 10.5 AU grow too quickly and migrate 
too fast to stay outside of 0.1 AU, indicating that planets that 
form in this part of the disc will end up as hot gas giants. When 
these planets reach r < 0.1 AU, the total evolution of the planet 
is stopped, explaining the kink in the gas mass of planets inside 
the grey area in Fig. [2 

Planets that form in the inner regions of the disc ro < 7 AU 
only grow a very small planetary core even though enough peb¬ 
bles are available, because the pebble isolation mass is low. The 
low core mass then prevents a fast accretion of the envelope, 
so that a long time is needed for the contraction of the enve¬ 
lope. This prevents planets that formed with ro < 5 AU from ac¬ 
creting a massive gaseous envelope, and these planets stay with 
Me > Menv, indicating that these planets are ice giants. These 
planets are also caught in a region of outward migration, letting 
them stay a few AU from the central star (Fig. [T]i- The planets 
that are within 5 AU and 7 AU have a slightly more massive 
core, indicating a slightly faster contraction phase of the enve¬ 
lope, so these planets can rapidly accrete gas and thus form gas 
giants that have then outgrown the region of zero migration and 
migrate towards the star. 

Our model predicts the formation of different types of planets 
by their initial formation location tq. Additionally, our model 
predicts that gas giant planets that have rf > 1 AU form in the 
outer regions of the disc (ro > 13 AU) an d thus do not form i n 
situ. This is a big contrast to the study of ICossou et al.l (l2014ll . 
where the cores of giant planets are built locally at a few AU 
by the accretion of planetesimals and planetary embryos. When 
the planetary cores then become massive enough, they can be 
caught in a region of outward migr ation (Fig.fll) and then accrete 
gas to form a gas giant. However, ICossou et al.l (l2014l) did not 
include gas accretion and therefore did not observe the inward 


migration of giant planets in type-II migration, which can bring 
them very close to the host star, especially when they form just 
a few AU away from it. In our study we overcome this problem 
because pebble accretion is very efficient in the outer disc and 
thus allows the formation of planets at large orbital distances, 
which allows them to stay far away from the host star even after 
their inward migration. 

3.3. Variation in the initial time and position 

We now expand Fig. [3] in the dimension of initial time tg when 
we place the planetary seed inside the disc. We start planets in 
our disc from to — 100 kyr up to the end of the disc’s lifetime 
of to = 3 Myr, as well as with 3 AU< ro < 50 AU. Figure |4] 
presents the final planetary mass as a function of initial radius ro 
and initial time to. Each point in this figure corresponds to the 
final mass of a growth simulation as in Fig. |2] All line cuts for 
a fixed initial time to with all 3 AU< ro < 50 AU correspond to 
a plot similar to Fig. |3] In particular. Fig. [3 represents a cut at 
fo = 2 Myr of Fig. |2 

In Fig. |4] the region ro and formation time to in combination 
with the final orbital distance rf can be interpreted as a map for 
the formation of different types of planets. Within rf <0.1 AU, 
we find hot gas giants. A late formation time to prolongs the 
formation time of the core, because fewer pebbles are available, 
and reduces the mass of the core at pebble isolation, because ///r 
drops in time. This means that the planet spends more time con¬ 
tracting its envelope, meaning that it can not accrete as much gas. 
Moreover, the late formation time reduces the time the planet mi¬ 
grates inwards in the disc, letting it stay farther out in the disc. 

In the band of 0.1 au< rf < 15.0 AU, we find warm and cold 
gas giant planets. The formation of these planets can also occur 
very easily until to ~ 2 Myr, because the pebble accretion rate is 
high even when the planets form at a large ro. This formation can 
then compensate for inward type-I migration through the first 
phases of gas accretion before the planets open up a gap in the 
disc and transitions into slow type-II migration. It also provides 
enough time to accrete a massive gaseous envelope after the long 
contraction phase. 

The planets around fo ~ 1 Myr and ro ~ 10 AU start in the 
flaring part of the disc, meaning that they have a have higher 
pebble isolation mass than the planets with smaller ro, so they 
can accrete their envelope faster and become more massive. At 
the same time, these planets are too massive to spend time in 
the region of outward migration, because they are already too 
massive when they reach this region (Fig. [T]i, which results in 
a continuous inward migration of these planets. They therefore 
reach 0.1 AU before they become more massive than ~ 300 
Earth masses. 

Eorming at a later time or an even larger ro results in smaller 
gas giants that are more in the mass regime of Saturn. Saturn- 
mass planets form in a distinct band below the white line in 
Eig. m where everything below the white line has Me < Menv, 
and runaway gas accretion has set in. All planets between the 
blue and the white lines in Eig. |4] have reached their pebble 
isolation mass and started to contract their envelope, but have 
Me > Menv In the inner disc these planets are only a few Earth 
masses and qualify as ice giants (light blue background colour in 
Fig. |4]between the blue and white lines), while in the outer disc 
the planets become more massive (Mp > lOMp) and are larger 
ice giants. The implications of our model regarding the forma¬ 
tion of the ice giants in our own solar system are discussed in 
more detail in section |5] 


Article number, page 9 of 123 1 


























A&A proofs: manuscript no. Planetgrowth 


3.0x10® 


2.5x10® 


2.0x10® 


1.5x10® 


1.0x10® 


0.5x10® 


5 10 15 20 25 30 35 40 45 50 

ro [AU] 

Fig. 4. Final masses of planets (total mass Mp = + M^nv) as a function of formation distance ro and formation time to in the disc. Planets that 

are below the dark blue line have reached pebble isolation mass and can accrete gas. All planets that are below the white line have < Me„v, 
indicating that they have undergone runaway gas accretion. The yellow line marks the pebble production line; planets below the yellow line cannot 
have formed by pebble accretion, because the pebbles have not yet formed at the insertion time of the planet. The black lines indicate the final 
orbital distance Pf of the planet. Each point in (ro, to) corresponds to the final mass of one individual evolution track as shown in Fig.|2] The white 
crosses indicate the initial orbital positions and times for the planetary evolution tracks shown in Fig.|^ 



In the late stages of the disc (t > 2 Myr), the region of out¬ 
ward migration is located at ~ 3 AU and can hold planets of up 
to a few Earth masses (Fig. [T]i. This means that ice giants stop 
their inward migration there, explaining the pile-up of ice giants 
in this region of parameter space. 

We want to note that Fig. |4] does not give a percentage of 
what kind of planets should exist around other stars. Figure |4] 
instead shows what kind of planet would form if an initial seed 
planetesimal with Mtrans (eq. E) was placed in the disc at ro and 
fo for a disc that lives, in total 3 Myr, around a solar type star. 
We can learn from Fig. |4]that gas giants are more abundant if an 
early formation scenario is invoked and that hot gas giants have 
to form early and fairly close to the central star, while cold gas 
giants can form at later times and further out in the disc. In situ 
formation of gas giants is not possible. For later formation times 
(fo > 2 Myr), the final mass of the planet becomes lower, be¬ 
cause the isolation mass becomes lower as Hjr drops, resulting 
in smaller planetary cores. A smaller planetary core then pro¬ 
longs the contraction time of the gaseous envelope, allowing the 
formation of planets that have Me > Menv at 2 Myr, because the 
planets have a shorter total evolution time. Considering that the 
most common observed exoplanets are small (Mp < lOMp), a 
later formation time for planets is favoured by our planet forma¬ 
tion scenario. 

These results emphasise the importance of an evolving disc 
structure that is not a simple power law. We discuss planet for¬ 
mation via pebble accretion in the MMSN in AnnendixICl where 


we show that the formation of different planetary types is dra¬ 
matically different when adopting too simplistic a disc model. 

3.4. Influence of pebble surface density 

The metallicity in pebbles Z has a strong influence on the out- 
come of planetary systems (iFambrechts & Johansenll2014t) . A 
higher metallicity in pebbles will allow for faster growth of the 
planetary core, resulting in a larger core. However, the pebble 
isolation mass (eq. reduces in time as the disc evolves and 
planets form, keeping the final core mass from being twice as 
high for discs with twice the metallicity in pebbles. Fower metal¬ 
licity in pebbles slows down the growth of the core of the planet, 
resulting in lower core mass, because the isolation mass drops in 
time owing to ///r decreasing in the longer time that is needed 
to build the core. This change in the mass of the core will then 
influence the gas accretion onto the planet, because a larger plan¬ 
etary core can contract its envelope in a shorter time than a 
smaller planetary core. The mass of the cores of these planets 
in a Z = 1.5% disc is ~ 20 - 30Me, which is about a factor of 
three to four too low to explain the amount in heavy elements 
for Corot-13 b, 14 b, 17 b, and 23 b, which have solid cores of 
around IOOMe (iMoutou et al.ll201^ . 

In Fig.|5]we display the final planetary masses for Z = 0.5% 
(top, half the nominal metallicity in pebbles) and for Z = 1.5% 
(bottom, 3/2 times the nominal metallicity in pebbles). In gen- 
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eral, we still find the same planetary classes in the ro - fo plane, 
but its distribution is different. 

In the low-metallicity case, the maximum to for forming a 
gas giant planet is reduced compared to the Z = 1.0% case, so 
that for to > 2 Myr, no gas giants can form any more. The region 
of parameter space that results in hot gas giants (rp <0.1 AU) 
extends to larger ro for a given to- This is caused by the fact that 
the planets in a disc with Z = 0.5% have smaller planetary cores, 
which results in a longer contraction time for the envelope. The 
planet therefore spends a longer time in the fast inward type-I 
migration regime, before it grows big enough via gas accretion to 
open a gap in the disc, which in total results in a smaller rp. This 
is another indication that 2 Myr is a reasonable seed formation 
time, because otherwise we observe more hot gas giants around 
lower metallicity stars. 

This effect not only influences the parameter space that har¬ 
bours hot gas giants, but it also influences all the parameter space 
where gas giants are formed. As a result, the region of parame¬ 
ter space with gas giants with 0.1 AU < rf < 5 AU is much 
smaller in the Z = 0.5% case than in the Z = 1.0% case, in¬ 
dicating that the gas giants should be rarer in discs with lower 
metallicity. Likewise, the region of parameter space for forming 
ice giants is slight ly enlarged in the Z = 0.5% case, consistent 
with observations (iBuchhave et alJiToidll . 

In the case with higher metallicity, Z = 1.5%, the parameter 
space for gas giants with rf < 0.1 AU is a bit smaller than in 
the Z = 1.0% case. Higher metallicity helps to keep planets out¬ 


side of 0.1 AU because the formation time of the core is shorter 
and the core is larger, also making the phase of envelope con¬ 
traction shorter, which in turn allows rapid gas accretion at an 
earlier stage. This then results in a planet being able to open a 
gap soon and migrate slower in type-II migration speed, allow¬ 
ing for larger rf in the end. Along with that, the more massive 
planets slow down their type-II migration rate because of the 
feedback with the disc (Ea.l22l). However, for later initial times 
fo > 1 Myr, this effect does not matter that much compared to 
the Z = 1.0% disc, because the planets do not reach such high 
masses that the reduction of the type-II migration speed plays 
that much of a role. Another aspect is that the formation of gas 
giants is now possible up to to = 2.5 Myr at nearly all initial 
orbital distances ro, except for the inner regions of the disc. 

On the other hand, the high metallicity in pebbles reduces 
the ro - to parameter space for the formation of ice giants slightly 
compared to simulations with low Z. A higher pebble accretion 
rate allows for efficient formation of planetary cores, even at late 
stages (high to) of the disc evolution, resulting in a more massive 
core than in the Z = 1.0% case, in turn allowing for a faster 
contraction of the envelope and thus Me < Mgnv The higher 
metallicity in pebbles then also allows for the formation of ice 
giants far out at late times in the disc, which was not that easy 
compared to the Z = 1.0% case. 

Generally, the higher metallicity significantly broadens the 
parameter space that allows for forming gas giants. A high peb¬ 
ble accretion rate is required to form massive gas giants (Mp > 
Mjup) that stay far outside in the disc (rf > 10 AU), which do not 
exist in the lower metallicity simulation. 

Observations also indicate that stars wi th higher metallic- 
ity host more giant plan ets in close orbits (ISantos et al.ll2004t 
iFischer & Valentil l2005h . which is confirmed by our results, 
where the formation of gas giants is possible in a wider range in 
the ro - to parameter space for high Z. If planet formation starts 
late (large fo), higher m etallicity helps to f orm giant planets at 
late sta ges. Additionallv. lMarcv et al.l(l2014l) and iBuchhave et all 
(12014ll find that smaller planets {Rp < 4f?E) are slightly more 
common around stars with solar metallicity, which is reproduced 
by our lower metallicities simulations, where the formation of 
small planets in the ro - fo parameter space is enhanced (top in 
Fig .|5]l compared to the high Z simulations. 

Fven with changing metallicities in pebbles that allow for 
faster or slower growth of the core, gas giants are not able to form 
in situ, because of their strong migration. On the other hand, the 
formation of ice giants results in much less migration through 
the disc, so these planets can form in situ, but they have to form 
at late stages, because otherwise they would continue to grow to 
become gas giants. 

3 . 5 . Amount of micrometre-sized dust in the disc 

The thermodynamic structure of the disc is determined through 
micrometre-sized dust grains. However, during the lifetime of 
the disc, the amount of micrometre-sized dust can change, for 
example because of grain growth. Larger grains do not con¬ 
tribute to the opacity, so that the smaller number of dust grains 
in the disc results i n a colder disc, because cooling is increased 
(iBitsch et al.ll20L5h . This results in a smaller aspect ratio of the 
disc, which reduces the pebble scale height (eq.lH and allows an 
earlier transition from 3D to 2D pebble accretion. Additionally, 
a smaller aspect ratio results in a larger Tf, making the transi¬ 
tion to 2D pebble accretion even smaller. Therefore planets in 
discs with a smaller amount of micrometre-sized dust grains can 
grow faster. However, a smaller aspect ratio of the disc reduces 
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Fig. 6. Final masses of planets (Mp = Me + Menv) as a function of initial 
radius tq and initial time fo in the disc. The plot features a metallicity 
of Z = 1.0% in pebbles and an amount of Zdust = 0.1% in micrometer 
sized dust grains, which is a factor of 5 smaller than in Fig.|4l The dif¬ 
ferent lines inside the plots have the same meaning as in Fig.|4l A lower 
amount of dust grains results in a colder disc, which reduces the pebble 
isolation mass in the inner disc and therefore increases the parameter 
space that allows for the formation of ice giants compared to Fig. [4] In 
the outer disc, formation of larger bodies is easier compared to Fig. [4] 
because of the larger pebble size and smaller pebble scale height /fpeb 
in the outer disc. 



the pebble isolation mass (eq. fTril) resulting the outcome of the 
systems. 

We now discuss simulations with an amount of micrometre¬ 
sized dust of Zdust = 0.1%, which is five times smaller than in our 
nominal model (see Fig.|4]i. All other parameters are the same as 
in our nominal model. Figure |6] shows a much larger parameter 
space in the inner parts of the disc that harbours ice giants than 
in the Zdust = 0.5% model. This is caused by the smaller aspect 
ratio in the inner di sc, which resulted from the increased cooling 
dBitsch et alJl2015h . This smaller aspect ratio results in a lower 
pebble isolation mass causing a longer contraction time of the 
gaseous envelope, which hinders the planets from reaching the 
runaway gas accretion stage. 

Changing the disc structure also influences the migration of 
embedded planets. In Fig. 0] the region harbouring ice giants at 
~ 3 AU is caused by a region of outward migration, which can 
contain planets of up to ~ 8 Mg (Fig.[T]i. However, in a disc with 
decreased Zdust the region of outward migration can only contain 
planets of up to ~ 3.5 Mg. The size of the cores in this region 
is only ~ 2.5 Me owing to the low pebble isolation mass. The 
cores therefore need to contract a gaseous envelope of ~ 2.5 Me 
before they can undergo rapid gas accretion, which takes a very 
long time. This means that the planets can outgrow the region of 
outward migration, because they become too massive, while they 
still undergo a contraction phase of the envelope. They therefore 
migrate to the inner system as ice giants, making the formation 
of ice giants in the inner system (rf <1.0 AU) possible, in con¬ 
trast to the Zdust = 0.5% model. This is crucial because most 
exoplanets that have been observe d have a few Earth m asses and 
orbit very close to the central star (iFressin et al.ll2013h . This area 
of parameter space can easily be populated in the case of low 
Zdust, which corresponds to cold protoplanetary discs. 

In the outer parts of the disc, planet growth seems to be more 
efficient than in the Zdust = 0.5% model for two main reasons. 
The increased cooling results in a smaller aspect ratio in the outer 
parts of the disc, which increases the size of the pebbles (Eq.O. 


Additionally, a larger pebble size reduces the scale height of the 
pebbles //peb. Both effects reduce the planetary mass needed to 
transition into the 2D pebble accretion branch (ea.fTOli. resulting 
in a faster growth rate compared to Zdust = 0.5%. This allows 
the efficient formation of gas giants in the outer disc, which then 
migrate into the inner disc. However, these planets will then only 
have a smaller planetary core than in the Zdust = 0.5% model, be¬ 
cause the pebble isolation mass reduces as well with decreasing 
aspect ratio (eq. \TB. 

A larger amount of micrometre-sized dust does not influ¬ 
ence the outcome of our simulations significantly, because the 
changes in th e disc structure are not as pronounced as for lower 
metallicities (iBitsch et alJl2015l) . 


4. Planet formation via planetesimal accretion 

In classic models of planet formation, planets grow via the ac¬ 
cretion of planetesimals. The isolation mass for planetesimal ac- 
cretion is different than when accreting pe bbles. It is given by 
(iKokubo & Idall200^ [Raymond et al.ll2014h 


Mi 


iso,pla 


■^^\I0 /?h) i 10 ) llAu) \Mo 


-0.5 


Me, 

(25) 


where b is the orbital separation of the growing embryos, which 
we set to 10f?H- Here, Zpia is the surface density in planetesimals, 
and Vpia the negative gradient of the surface density in planetes¬ 
imals. The accretion rate by plan etesimal accretion is cha nged 
compared to pebble accretion. In iLambrechts et al.l (l20I4l) . the 
accretion rate for planetesimal is given by 


Afc.plan “ 'Filifc.peb — 'FridOpjZpeb , 


(26) 


where is a reduction factor to the normal pebble accretion rate. 
It is given by 


= 3 X 10 


UOAU 


(27) 


This follows directly from the assumption that 

the planetesimal veloci ty dispersion is equal _ to 

the Hill speed ( Dodso n-Robinsonet al.l 120091 : 
iDodson-Robinson & Bodenheimei 201 Oh . hh = D^h, and 
gravitational focusing occurs from a radius (rcrH)*^^rH, which is 
smaller than the planetesimal scale height F/pia = va/^ - ra- 
The general assumption of simulations with planetesimal ac¬ 
cretion is that all solids in the disc are turned into planetesimals 
at the start of the simulation. We do the same here and point out 
that the surface density in planetesimals is Epeb = ZEg, where 
we set Z = 8.0%, which is eight times higher than for the accre¬ 
tion with pebbles, and about eight times higher than the assumed 
metallicity of 1% in the MMSN. With Z = 1.0% planetesimal 
accretion is too slow, and no planets with Mp > IMe are formed 
at all. Along with that, we keep the surface density of planetes¬ 
imals constant through the evolution time of the disc, because 
planetesimals are safe from drifting through the disc by gas drag. 
This also means that we keep the initial gradients in the surface 
density of planetesimals that follow the initial gradients of the 
gas surface density Zg. When the simulation starts at a later ini¬ 
tial time fo, we use a planetesimal surface density and gradient 
that follows Zg at that time fo. When the planets reach their plan¬ 
etesimal isolation mass, they first have to undergo a contraction 
of the envelope before runaway gas accretion can start. However, 
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Fig. 7. Final masses of planets forming by planetesimal accretion as 
a function of initial radius ro and initial time fo in the disc for Z = 
8.0% in planetesimals. The green line is the pebble production line, 
where we assume that a production of pebbles is needed to first form 
the planetesimals, so every point beneath the yellow line should not be 
taken into account. The red line now indicates the planetesimal isolation 
mass fea. l26b . so everything below it has reached planetesimal isolation 
mass. Everything above the white line has The top plot 

shows the full ro-to parameter space, while the bottom plot is zoom into 
the inner rg-to parameter space to enhance the interesting parts of the 
diagram. 


lation none of the planets reach a mass that is comparable to 
Jupiter’s. Additionally, these simulations fail to produce planets 
with Mp > IOMe outside of 8 AU. 

With just the accretion of planetesimals it is very hard to 
form the cores of gas giants that stay out at ~ 5 AU, as does 
Jupiter in our own solar system. However, if multiple planetary 
embryos were present in the disc, these could collide and form 
bigger objects an d eventually the cor es of giant planets at these 
orbital distances dCossou et al.ll2014t) . With single embryos that 
accrete the planetesimals, effective growth to reach the stages of 
giant planets is not possible at all. 

The reason we do not produce giant planets at a few AU in 
contrast to population synthesis models lies in their simpler disc 
model. The population synthesis models use a steeper gradient 
in planetesimal surface density , allowing for mo re planetesimals 
in the inner parts of the disc (llda & Linll2008h . A higher den¬ 
sity of planetesimals then results in a faster growth rate of the 
cores, which also allows for a faster contraction of the gaseous 
envelope, giving the planet more time to reach the runaway gas 
accretion stage. Still, formation of gas giants in orbits beyond 
5 AU would not be possible, even under these generous condi¬ 
tions. 


5. Formation of the solar system 

We now focus on the formation of the giant planets in our solar 
system, where we follow two different approaches. We first want 
to repr oduce the planetary c onfiguration at the start of the Nice 
model (iTsiganis et al.ll200^ . and in a second attempt we want to 
reproduce the giant p lanet configuration at the beginning of the 
Grand Tack scenario (IWalsh et al.ll201 Ih . where we end our sim¬ 
ulations just before the two gas giant planets start their outward 
migration in resonance. We stick here to our usual assumption 
of a metallicity of 1.0% in pebbles that can be accreted onto the 
planet. 

When multiple planets form in the disc, the outermost pebble 
accreting planet will reduce the flux of pebbles seen by the inner 
planets by the amount that it accretes. Howeyer, this reduction 
of the pebble accreti on stream is not signific ant when just con¬ 
sidering four bodies dUambrechts et al.ll2014l) . 


since the time scale of forming cores with planetesimals is yery 
long, we extend the disc’s lifetime to 5 Myr to extend the time 
for gas accretion. This is also because a disc lifetime of 3 Myr 
did not produce any gas giants whatsoeyer. 

In Fig. |7] we display the final planetary masses for simula¬ 
tions where the planetary cores grow by the accretion of plan¬ 
etesimals with a metallicity of planetesimals of Z = 8.0%. The 
growth of planets, eyen with eight times higher metallicity is 
yery slow, and only those planets that start to form early (small 
fo) and in the inner regions of the disc (small ro) reach planetesi¬ 
mal isolation mass. Howeyer, the planets that reach planetesimal 
isolation mass are yery small and take a yery long time to con¬ 
tract their enyelope, which means they will migrate for a yery 
long time in type-I migration. But the planets are indeed small 
enough, and in the inner regions of the disc these planets can 
get trapped in the zero-migration zones (Fig. [T]) allowing them 
to stay at a few AU from their host star before runaway gas ac¬ 
cretion sets in. The planets then outgrow the region of outward 
migration and moye towards the inner disc. In the end most of 
these planets are not massiye enough to open up a gap in the 
disc, and they migrate inwards yery fast. In fact, in this simu- 


5.1. Nice model 

The Nice model aims to explain the bombardment history of the 
inne r solar system, rough ly 1 Gyr after the formation of the sys¬ 
tem (iTsiganis et al.lUOOa) . In this model, the four giant planets 
in the system start in compact orbits after the gas disc dispersed, 
where Jupiter and Saturn are in a 3;2 resonance, and outside the 
ice giants is a belt of planetesimals containing roughly 30 Earth 
masses. The system then becomes unstable, because of the con¬ 
stant scattering of planetesimals on the planets. This instabil¬ 
ity leads to an inward motion of Jupiter, while the other planets 
moye outwards. In ab out 50% of the simu lations, Neptune and 
Uranus switch places (ITsiganis et al.ll20d^ . 

We now aim to reproduce the masses and planetary orbits of 
the four giant planets in the solar system using our disc eyolution 
and pebble accretion model. We keep the total disc lifetime at 
3 Myr, but also allow different formation times for the planetary 
seeds of the giant planets. The results are shown in Fig.[3 which 
shows the final eyolution of the system at fo = 3 My. The starting 
times of the planetary seeds differ between the planets. The final 
planetary masses and orbits needed for the starting configuration 
of the Nice model are reproduced quite well. 
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Fig. 8. Evolution of the total planetary masses as a function of orbital 
distance for a proto-solar system with the giant planets ending up in a 
configuration similar to the initial conditions of the Nice model. Jupiter 
starts at to = 1.65 Myr, Saturn at to = 1.77 Myr, and Neptune and 
Uranus both at to = 1.81 Myr. The big black circles label a total disc 
evolution time of 3.0 Myr, while the small black dots indicate 2.0, 2.2, 
2.4, 2.6, and 2.8 Myr. We display Neptune inside of Uranus here, sinc e 
they often switch places during the Nice model dTsiganis et alJl2Q()^ . 
We used - 0.5% and a metallicity of pebbles of Z = 1.0% as in 
our nominal disc model. 


The formation of a gas giant is easy in the framework of the 
evolving disc and pebble accretion. However, the formation of 
the two gas giants of the solar system is more somewhat trickier, 
because in the final configuration, they differ by a factor of 3 in 
mass, but are only 3 AU apart from each other. The reproduction 
of this exact conhguration is challenging because either Jupiter 
is slightly too small or Saturn is too big. This is a consequence of 
the interplay between pebble accretion, planet migration, gas en¬ 
velope contraction, and runaway gas accretion. In the simulation 
shown here, our Jupiter analogue is only ~ 230 Earth masses. 
This difference compared to the real Jupiter can be caused by the 
slight uncertainties in our model, for example in the contraction 
of the gaseous envelope (see Appendix I bTi. Even though there is 
a time difference between the formation of Jupiter and Saturn, 
their orbits also start very close to each other, which might result 
in some minor interactions between the planets that is not taken 
into account here. Nevertheless, the reproduction of Jupiter’s and 
Saturn’s masses and orbital distances is remarkably good. 

The formation of the ice giants that are in the mass range of 
Uranus and Neptune with a very low gas content is difficult be¬ 
cause a larger core can attract gas more easily (eq. [nil, so the 
planet can very easily grow to become a gas giant. This limits 
the parameter space in ro-fo that hosts ice giants with low gas 
content in the mass range of Uranus and Neptune at hnal or- 
bital distances of 15 -20 AU (Eig.lH). However, in a recent study, 
llzidoro et alJ (1201 5h have found that the formation of Uranus and 
Neptune by a series of giant impacts of planetary embryos of 
a few Earth masses outside of Saturns’ orbit is possible. These 
planetary embryos can easily be formed by seeds growing with 
pebble accretion. When taking also the possibility of giant im¬ 
pacts into account, the formation of ice giants like in our so¬ 
lar system seems very likely, especially at late formation times. 
Eull-grown ice giants and the planetary embryos can form in or¬ 
bital distances of 20 AU to 30 AU at 2 Myr (Eig.|3]). The region 
where ice giants can form in the ro-fo parameter space is also 
indicated in Eig. [TO] It is a vast region in the parameter space. 


We also checked that a change in Zdust influences our results 
when reproducing the initial conhgurations of the Nice model, 
but found no such dependency. Reproducing the exact conhgu¬ 
ration of the solar system is very sensitive to the initial orbital 
distance and formation times of the planetary seeds, since those 
parameters determine the hnal planetary mass and orbital dis¬ 
tance. The general outcome of having two gas giants, and outside 
of them, two ice giants can be reproduced very easily. 

The orbital conhguration of the solar system in the present 
day is different from the initial conditions of the Nice model. 
The ice giants are now much farther away from the Sun than in 
the initial conhguration of the Nice model. In particular, Neptune 
is located 30 AU from the Sun. Eorming a Neptune-sized body 
with rf = 30 AU in our planet growth scheme is possible, but the 
planetary seed would then have to start deep within the Kuiper 
belt at ~ 45 AU and with an early formation time to ~ 670 kyr. 
When the planetary seed reaches 40 AU, it has already grown 
to several Earth masses, which would have disrupted the Kuiper 
belt. This supports the Nice model concept of forming the giant 
planets in a close resonant conhguration, because it leaves the 
Kuiper belt untouched. 

5.2. Grand Tack scenario 

The Grand Tack scenario (I Walsh et al.ll201 ih describes a sce¬ 
nario where Jupiter and Saturn migrate into the inner solar sys¬ 
tem (Jupiter down to ~ 1.5 AU) and then migrate outwards in 
resonance again. This outward migration o f gap-opening planets 
in reso nance was originally discovered bv iMasset & Snellgrovel 
(1200Ih . but is applied to the solar system here. The appealing 
effect of the Grand Tack scenario is that the masses and orbital 
distances of the terrestrial planets, especially Earth and Mars, 
can easily be repro duced, as can the features of the asteroid belt 
(I Walsh et al.ll201 ih . Here we now want to reproduce the masses 
and orbital distance of Jupiter and Saturn before both planets be¬ 
gin their outward migration in resonance. 

Because we only model single planets in the disc, we try 
to arrive at a conhguration that would allow for the outward 
motion of the giant planets in the gas disc in resonance. Addi¬ 
tionally, because the outward migration happens in a gas disc, 
the planets will still accrete gas during this outward motion. We 
therefore aim to have a Jupiter analogue with ~ 200 Me and a 
Saturn analogue with ~ 70 Me. This is much smaller than their 
current mass, but during their outward migration in the gas disc, 
Jupiter and Saturn can still accrete gas and reach their hnal mass. 
We also aim for stranding the planets roughly in a mutual reso¬ 
nance, where we put Jupiter at ~ 2 AU. We try to put Saturn’s 
orbit somewhere between the 2:1 and 3:2 resonance with Jupiter, 
because b oth resonances allow outward migration in discs with 
small Hj r (iPierens et al .12014 . as we have here. We also assume 
that when the disc reaches 3 Myr of lifetime, the disc will still 
live long enough to allow the giant planets to migrate outwards 
and accrete the rest of their mass. 

Eigure 0 shows the ro-to parameter space of interest for the 
formation of Jupiter and Saturn with the above-mentioned char¬ 
acteristics in a disc with Zdust = 0.1%. The initial ro of Jupiter 
and Saturn is quite close to each other, but they are separated 
in to by nearly 200 kyr. During this time, Jupiter’s core grows 
to a few tenths of an Earth mass and migrates only a minimum 
distance through the disc. However, before the mass of Jupiter 
becomes massive enough to disturb the orbit of Saturn’s seeds, 
it migrate a few AU from its initial position. _ 

In the Grand Tack scenario proposed by IWalsh et"^ (1201 Ih 
Saturn migrates much faster than Jupiter and catches it. This is 
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Fig. 9. Configuration of our assumptions of the Grand Tack scenario 
in a disc with = 0.1%. The black lines indicate the orbital dis¬ 
tances of 2.0 and 3.5 AU, which correspond to the 2 : 1 mean motion 
resonance when the outward migration in the Grand Tack scenario can 
start. The white lines correspond to the final planetary masses of 200 
Me for Jupiter and 70 Me for Saturn. Jupiter’s rg and tg are therefore 
constrained by the crossing between the white 200 Me and the black 
2 AU lines, while Saturn is constrained by the white 70 Me line and the 
black 3.5 AU line. Jupiter forms earlier (smaller tg) and closer to the 
Sun (smaller rg) than Saturn, showing that the conditions of the Grand 
Tack scenario can be met. 



explained by type-III migration. However, our Jupiter and Sat¬ 
urn seeds are both placed in a region of parameter space that is 
subject to type-III migration (see appendix|0. This indicates that 
Jupiter must stop its inward type-III migration and transition into 
type-II migration, so that Saturn catches up. However, this pro¬ 
cess is not investigated in full. Likewise, the outward migration 
of Jupiter and Saturn in resonanc e is thought to end when th e 
aspect ratio of the disc increases dMasset & Snellgrovell200ll) . 
However discs with low numbers of dust grains at the late evo¬ 
lutionary stages have low aspect ratios t o greater distances than 
Saturns current orbit (iBitsch et al.ll2015h . 

However, the formation of Jupiter and Saturn for the con¬ 
figuration of the Grand Tack scenario is only possible in discs 
with a small number of dust particles. We did not find not find 
any allowed parameter-combination in ro - to space for discs 
with large amounts of micrometre-sized dust (Zdust > 0.3%) that 
would allow for the starting configuration of the Grand Tack sce¬ 
nario within our assumptions above. Given the uncertainties in 
our model (e.g. gas envelope contraction), more investigations 
are certainly needed to reproduce the initial configuration of the 
giant planets in the Grand Tack scenario. 


6. Discussion 

6.1. Evolution and lifetime of the disc 


migrate significantly in the disc, because migration is mostly in¬ 
wards at that stage of the disc’s lifetime. Combined with the long 
time scale of 5 Myr, planets can only stay outside of 0.1 AU if 
they for m late (after ~ 2.5 Myr) and far outside in the disc. In the 
study bv iMamaiekI (l2009h only ~ 30% of all stars are found to 
still harbour protoplanetary discs after 2.5 Myr, so our disc life¬ 
time of 3 Myr might be even on the longer side. Different disc 
lifetimes are discussed in more detail in AppendixlAl 

The evolution and formation of planets inside the disc de¬ 
pends crucially on the disc structure and its evolution. If the disc 
can sustain higher M rates for a longer time, planet formation is 
significantly different, because the accretion rates will be higher 
and the disc might entertain regions of outward migration for 
higher planetary masses for a longer time, allowing more planets 
to stay farther out in the disc. Improved observational constraints 
on the lifetimes of protoplanetary discs would therefore be very 
useful for understanding planet formation better. 


6.2. Inner disc 

We stop the inward migration when the planet reaches rf < 
0.1 AU, which is the inner edge of the disc. We then stop the 
whole evolution of the planet at that point, because tidal effects 
so close to the star can become important, which we do not con¬ 
sider here. 

In additional, planetary growth that close to the central star 
might work differently; for example, the pebbles that can be ac¬ 
creted from planets that are at rp < Hce are significantly smaller 
than for rp > Hce. Silicate particles w ill bounce off each other in¬ 
stead of growing (IZsom et al.ll20K)l) . keeping the particles small. 
We also stop the growth of the planet, because gas accretion in 
the inner hot parts of the disc might work completely differently 
from gas accretion in the outer colder parts of the disc, because 
the surrounding gas temperature is higher, making cooling less 
efficient, thereby prolonging the contraction phase of the enve¬ 
lope. The final mass of the planets inside the parameter space ro 
and to, which allows for rf < 0.1, is therefore not 100% accurate, 
but orbital evolution is taken into account self-consistently until 
the planets reach 0.1 AU. 

An additional challenge emerges when looking at results 
from N-body simulations, where the inner edge of the disc acts 
as a halting zone for inward migration. In ICossou et akl (12014h 
the inner edge acts as small zone of outward migration that stops 
planets. These planets can then act as a buffer zone for other 
embryos migrating into the inner disc and allow pile-ups in res¬ 
onances. The planets trapped in resonance can grow more and 
become bigger planets farther out in the disc. Also, in our simu¬ 
lations, only gas giants migrate to sub-AU orbits, while smaller 
ice giants remain around their birth locations. 

6.3. Convergence zones 


The evolution of the disc is c oupled to the acc r etion rate, where 
we follow the observations of lHartmann et alJ (Il998h . We make 
the assumption that when the accretion rate is very low, M < 
2 X 10~^Mo/yr, photo evaporation clears the disc very rapidly 
([Alexander et al.ll2014t) . and we end our simulation at that stage. 
This results in a lifetime of 3 Myr. If the disc evolution contin¬ 
ues to M = lx lO^^Mo/yr, the lifetime would be 5 Myr. This 
additional lifetime of 2 Myr influences the results because in the 
late stages of the disc, the disc mass is low so the pebble ac¬ 
cretion rate onto the planet is low as well. The planets therefore 
no longer grow efficiently, but at the same time the planets can 


In the inner regions of the disc, where Hjr drops with radius, 
zones of outward migration exist. In the early stages of the disc 
evolution, these zones of outward migration are several AU wide 
and can harbour planets from 5 to 40 Earth masses. However, as 
the disc accretes in time, the regions of outward migration shrink 
and can only contain planets with lower mass (see Fig.ITJ. At an 
age of ~ 1 Myr the regions of outward migration exist only in 
the inner few AU and for planets with Mp <10 Me. 

This implies that planets that form early in the disc will out¬ 
grow the regions of outward migration owing to gas accretion 
and migrate to the inner disc. In the early stages of the disc evolu- 
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tion, the region of outward migration allows the planets to grow 
to several tens of Earth masses before they migrate to the inner 
disc, but it does not prevent the planets from migrating to the 
inner disc. It just delays their inward motion. 

In the late stages of the disc evolution, in combination with 
the slow gas accretion due to the long contraction time of the 
envelope, however, the region of outward migration is crucial. 
Late-forming planets (fo > 2 Myr) in the inner disc have only 
a very small core, because the pebble isolation mass is low 
(~ 3Me). Therefore the contraction of the envelope takes a long 
time (eq. [nil, which hinders the planet from outgrowing the re¬ 
gion of outward migration and lets it stay at a few AU from the 
central star. These planets in the inner disc then classify as ice 
giants, because they are just a few Earth masses, their core mass 
is higher than their envelope mass, and they form at rp > rice. 

Additionally, the in situ formation of giant planets is not pos¬ 
sible when we take migration into account. Even planets that 
have a final orbital distance of a few AU have most likely not 
formed there, but migrated inwards from even farther out (see 
Eig.m. Smaller planets, on the other hand, can form in situ, but 
only i n orbits beyond 1 AU. 

In iLambrechts & .lohansenl (1201 4l) the planets form in an 
inside-out fashion, where planets closer to the star form before 
planets farther out, when the seed cores from around the same 
time in the disc. However, inside ten AU, this simple mechanism 
breaks down in our model. Time-evolving discs including con¬ 
vergence zones for migrating planets do not allow such a simple 
picture any more (see Eigs. 1 and|4]i. 

6.4. Formation of the solar system 

The formation of the exact configuration of our own solar sys¬ 
tem is relatively straightforward. The formation of the cores of 
the giant planets can eas ily be explained with pebble accretion 
(ILambrechts et al.ll2014l) . however reproducing the exact plan¬ 
etary masses and orbital distances requires some fine-tuning of 
the model. The two main problems are (i) that there are unknown 
opacities in the planetary atmosphere when the gaseous envelope 
contracts, which changes the time scale of gas envelope contrac¬ 
tion; and (ii) that we were just considering single planets in the 
disc, meaning that we do not c onsider giant i mpact s for the for¬ 
mation of the cores. A study bv llzidoro et al.l (1201 5h explains the 
formation of Uranus and Neptune through the accretion of plane¬ 
tary embryos of a few Earth masses, which can be formed easily 
by pebble accretion. Taking this growth mechanism into account 
makes the formation of the ice giants of our solar system even 
easier. 

Re producing the initi al conditions for the Grand Tack sce¬ 
nario dWalsh et al.l 1201 ih is more difficult because only discs 
with a low amount of dust particles (low Zdust) allow reproduc¬ 
tion of initial conditions that can lead to the Grand Tack sce¬ 
nario. The problem faced is the necessarily big mass difference 
between Jupiter and Saturn even though they are quite close in 
their orbital separation at the beginning of their combined out¬ 
ward migration. During their outward migration, Jupiter and Sat¬ 
urn can still accrete gas, which would allow them to grow, but 
making it very hard to constrain their initial mass at the start of 
their outward migration. This scenario should be investigated in 
more detail with simulations that feature multiple cores. 

The formation of the orbit al configuration need ed for the late 
instability in the Nice model dTsiganis et al.ll200^ can be repro¬ 
duced quite well within our simulations, where the amount of 
dust particles for the disc structure does not play an important 
role. The reproduction of the Nice model configuration is eas¬ 


ier because larger orbital distances between the giant planets are 
needed than in the Grand Tack scenario, which allows the mass 
ratio between the giant planets to be reproduced. This is much 
harder in the Grand Tack scenario (Eig. 01. The orbital distances 
of the giant planets are matched for our Nice model configura¬ 
tion, while only the mass of Jupiter is slightly off (see Eig. O. 
Our model is therefore very successful in reproducing the con¬ 
figuration of the giant planets in our solar system. 

6.5. Formation of Jupiter 

The atmosphere of Jupiter is enriched in noble gases compared 
to the solar values dOwen et al.lll999h . which led to the conclu¬ 
sion that the planete simals accreted by Jupiter formed in the cold 
regions of the disc dOwen et al.l 1 19991) . In this scenario the for¬ 
mation of Jupiter’s core starts at 5 AU, and planetesimals from 
the cold region of the solar system are accreted in order t o pro- 
vide the enrichments in noble gases. iGuillot & Huesol d2006l) 
proposed that the accretion of Jupiter’s noble gases occurs at a 
late stage of the disc evolution, probably during the times when 
photoevaporation of the disc already has started because the tem¬ 
perature is low and noble gases can condense out. 

However, our Jupiter analogue (Eig. Ell starts to form at 
20 AU and reaches pebble isolation mass when it has migrated 
down to ~ 17 AU. Since the formation starts at fo = 2 Myr, the 
surrounding disc temperature in that region is already very cold 
(T ^ 30K), which allows for the condensation of ices that en¬ 
gulf noble gases and which then would form the core of Jupiter. 
Moreover, the accretion of gas onto the core also happens in the 
cold regions of the disc, which could further explain the enrich¬ 
ment of noble gases in Jupiter’s atmosphere. Our model implies a 
late formation of Jupiter at a large distance and therefore would 
naturally explain the enrichment of noble gases in Jupiter’s at¬ 
mosphere. 

Our model therefore also predicts a noble gas enrichment 
in Saturns atmosphere, because Saturn forms in the cold parts 
of the disc. Unfortunately, the noble gas enrichment of Saturn’s 
atmosphere is unknown at the moment. 


6.6. Formation of hot gas giants 


The formation of hot gas giants close to the central star can 
happen t hrough two form a tion channels . Either by inward mi¬ 
gratio n (iLin et al.l 119961 : IWardI 119971: iMasset^Papaloizod 
20031) or by scattering ev ents (iRasio & Eordi 19961 : 
Marzari & Weidenschillingl l2002l) . Inward migration of gas 
giants through type-II or type-III is reasonable possibility in 
our formation scenario, because the migration rates are very 
strong (Eig. H. However, when comparing the simulations with 
Z = 1.0% and Z = 1.5% in pebbles, it seems that the amount 
of parameter space needed to form gas giants with ff < 0.1 AU 
is slightly smaller in the Z = 1.5% case. This effect is caused 
by the faster growth rate of planets, making the planets in 
the Z = 1.5% disc more massive than in the Z = 1.0% case, 
which then results in a reduced type-II migration rate because 
the planets become so big that they slow down their inward 
migration (eq. l2^ . However, this difference in ro-fo parameter 
space might be a bit deluded by type-III migration, which acts 
on more planets in the high Z case (see appendix |E]i, which in 
turn might lead to more planets reaching rf < 0.1 AU in the high 
Z case. 

Scattering events are also used to explain the distribution of 
hot gas giants, especially for hot gas giants with high eccentric- 
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ity dMarzari & Weidenschillingl2002HDavies et al.ll2014h . These 
scattering events can be caused by other giant planets far out in 
the system. This clearly favours the high Z case, because the for¬ 
mation of giant planets at large orbits (up to 30 AU) is possible 
easily until the disc has reached an age of ~ 2.0 Myr (Fig. |5]l, 
while for the lower Z cases it is hard to form giant planets even 
up to ~ 15 AU (Fig.m. The larger formation probability of gas 
giants at greater distances makes the scattering hypothesis more 
probable, because there can simply be more giant planets avail¬ 
able that can be scattered to the inner system to form a hot gas 
giant. 

The combination of those two possibilities (disc migration 
and scattering) supports observations that hot g as giants are 
found around host stars with higher metallicity (ISantos et al.l 
l2004t Irischer & Valentill2005[). Many hot gas gia nts have a large 
content of heavy elements JMoutouet aklboi^ . which can be 
explained by scattering events of planets that formed far out in 
the disc, because the core mass in the outer disc is higher owing 
to the higher pebble isolation mass. 

6.7. Planetesimal accretion 

The accr etion of planetesimals can be slower than the lifetime o f 
the disc (iPollack et al.l 1 995 lRafikovl2004tlLevison et'ani2010l) . 
even when considering an increased amount of solids inside the 
disc. We confirm these results (Fig.|7]i and show that pebble ac¬ 
cretion is much more efficient, even at the lower metallicity of 
pebbles than for planetesimals (Fig. Ell. In our nominal disc that 
has 0.5% of metallicity in micrometre-sized dust and 1.0% of 
metallicity in pebbles to form the planets, we are able to form a 
solar system analogue, which is not possible with planetesimals. 

Planetesimal accretion with a metallicity of 8.0% in plan¬ 
etesimals is only able to explain the formation of gas giants in 
the inner regions of the disc (rf < 1 AU), if along with the high 
metallicity, the lifetime is prolonged to at least 4 Myr. This is 
another limitation, because at that age p hotoevaporation m ight 
clear the disc efficiently. Observations bv iMamaiekI (l2009l) also 
indicate that only ~ 30% of all stars still harbour protoplanetary 
discs after an age of 2.5 Myr. In our simulations with planetesi¬ 
mal accretio n, we are not able to form a Jupiter planet at 5 AU, 
in contrast to lPollack et al.l (Il996h . because the isolation mass is 
lower there. This in turn leads to a longer contraction time for 
the envelope to prevent the accretion of gas, because our pro¬ 
toplanetary disc models are much colder than the MMSN. Ad¬ 
ditionally, we included the migration of planets, which leads to 
inward migration of planets, s o that it is no t possi ble to form a 
Jupiter analogue at 5 AU as in iPollack et akl (Il996l) . 

In addition, planetesimal accretion is not able to form the 
cores of the ice giants in the outer system. Even with an enhance¬ 
ment of the planetesimal density by a factor eight, planetesimal 
accretion alone can therefore not explain the formation of the 
solar system. 

6.8. Ice planets 

Our simulations predict a new class of planets, the so-called ice 
planets. They have a more massive core than envelope, but only 
reach a few Earth masses, so they do not qualify as ice giants. 
We also call them ice planets because they form outside the 
snow line (rp > rice), meaning they mostly consist of ices and 
some gas, so they should have a lower density than rocky plan¬ 
ets. These ice planets can be formed in the late stages of the disc 
evolution (fo > 2 Myr) in the outer regions of the disc. In these 


parts of the disc, the planets do not reach pebble isolation mass, 
because there is not enough material available at the late stages 
of the disc evolution, hindering growth. If formation of planets 
takes place generally at later stages in the disc evolution, and if 
the disc can still provide pebbles at these late stages, ice plan¬ 
ets should be very common and be detectable by observations 
at large distances around their host stars. These type of planets 
are not predic ted, however, and are not seen by population syn¬ 
thesis studies (iDittkrist et al.ll20f^ . because the formation with 
planetesimals takes too long, even for these low masses (Eig.|7]i. 

6.9. Population synthesis 

Population synthesis models attempt to reproduce the observed 
distributions of exoplanets by combining planetesimal accretion, 
gas accretion, and planet migration in protoplanetary discs. The 
first investigations of this nature sta rted nearly a decade ago 
(llda & Linll2004 lAlibert et ^l2004l) . However, to make these 
models work, some fundamental assumptions are made. 

Fi rstly, because the ac cretion of planetesimals is very ineffi¬ 
cient (iLevison et al.ll201^ . a much higher total amount of solids 
than the standard solar value has to be assume d in order to make 
the population synthesis simulations work. In llda & LinI (l2004h 
the total amount of planetesimals has to be at least five times the 
solar value in order to produce planets effectively. 

Secondly, the migration rates of pro t oplan ets inside the discs 
are found to be too high. llda & LinI (l2008l) find that type-I- 
migration needs to be reduced by factors of 10 -100 to reproduce 
the observed distributions of exoplanets. However, by taking 
outward migration due to the en tropy driven corotation torque 
(iPaardekooper & Mellemal[2006h into account, the reduction of 
type-I migration can be greatly decreased so that it might no t 
be needed any more (lAlibert et al.ll2013t IDittkrist et alfcofil) . 
However, in these si mulations ~ 60% of t he formed planets still 
move inside 0.1 AU (IDittkrist et al.ll20i4l) . showing that migra- 
tion is still a pro blem for planet population synthesis models 
(iBenz et al.ll2014t) . 

Thirdly, the evol ution of the disc i s either consider ed only via 
a simple power law (llda & Linl2004lida & Linll2008h or follows 
a viscous disc evolution, where the gradients in surface density 
and temperature do not change in value b ut only move inwards 
as the disc accretes (IDittkrist et al.ll20f^ . Keeping the gradient 
in surface density and temperature at the same value has impor¬ 
tant consequences for the torque acting on embedded planets in 
the disc, because these gradients determine the migration rates 
and direction of planets (see section 1274) 1. However, recent simu¬ 
lations of protoplanetary discs taking heating and cooling effects 
into account ha ve shown that these gradients change as the disc 
evolves in time (iBitsch et al.ll2015l:lBaillie et al.ll2015l) . 

Fourthly, the initial gas mass, metallicity, and lifetime of the 
disc and the initial starting position of the p lanetary seeds ar e 
randomized in population synthesis studies (iBenz et al.ll2014l) . 
This can allow for the formation of massive planets (Mp > 
IOOMe) in the outer regions of the disc (rp 10 AU), where 
the random selection of a very high disc mass allowed the for¬ 
mation of massive planets even in the outer disc. Disc column 
densities of more than ten times the M MSN are not found in ob - 
servations of real protoplanetary discs (I Williams & Ciezal201 ll) . 
However, even for such high disc masses, population synthesis 
models have problems forming planets at rp > 10 AU. 

After applying these fundamental assumptions to their mod¬ 
els, population synthesis models can reproduce the observed dis¬ 
tribution of exoplanets to some extent. However, they still face 
two challenges: i) a huge amount of the initial seed planets mi- 
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grate to the inner edge of the disc and ii) they fail to reproduce 
ice giant planets in the outer parts of the disc (rf > 8 AU) on a 
regular basis. These type of planets are only produced if the con¬ 
ditions for the surface density of solids and disc dispersal time 
are just right by coincidence. 

In the pebble accretion scenario we do not make any simpli¬ 
fications or reductions regarding the migration speed of planets, 
nor do we increase the amount of solids in our disc to unreason¬ 
ably high amounts. Our planets can form with the nominal solar 
metallicity, where we use 0.5% of solids in micrometre-sized 
dust grains and 1.0% of solids in the form of pebbles, which 
makes a total of 1.5% in solids. 

The amount of planets that reach the inner edge of the disc is 
greatly reduced when the planets start to grow late in the disc. 
Considering that only a few percent of the stars harbour hot 
gas giants, a late formation scenario of planetary seeds (to ~ 2 
Myr for disc lifetime of 3 Myr) with pebble accretion natu¬ 
rally explains this observation, because the available time un¬ 
til disc dispersion is reduced, making the formation of hot gas 
giants more difficult. The late formation of planets also results 
in the formation of ice giants in the inner disc (between 1 and 
5 AU). Those are among st the most detected planets in our 
galaxy (iFressin et alj|2013h . 

Pebble accretion is very efficient in the outer disc in the early 
stages (Fig. SI, but becomes less efficient as time evolves. How¬ 
ever, pebble accretion is much more efficient than planetesimal 
accretion, which allows the formation of ice planets in the outer 
disc. 

In contrast to population synthesis models, where the out¬ 
come (which type of planet) is determined by the initial condi¬ 
tions such as the mass and lifetime of the protoplanetary disc 
and the surface density of planetesimals, our pebble accretion 
approach shows that the outcome of planet formation can be 
very different depending on the different initial times at which 
the planetary seed is placed in the disc. 

We have also tested whe ther planets in our simulatio ns are 
subject to type-III migration dMasset & PaDaloizoull2003]) or are 
influenced by dynamical corotation torques (see Appendix|0. In 
our simulations, we find that these effects only play a role in the 
early stages of the disc evolution, when the disc is still massive. 
At the later stages of the disc evolution, these effects no longer 
play a significant role because the disc is not that massive any 
more. 

The effect of type-III migration and of the dynamical corota¬ 
tion torques are crucial for population synthesis models, where 
a higher disc mass is needed to form planets. However, popu¬ 
lation synthesis models do not take type-III migration into ac¬ 
count, which makes evolution of the semi-major axes in popula¬ 
tion synthesis models very doubtful, because type-III migration 
and dynamical corotation torques can lead to a faster inward mi¬ 
gration. This can result in even more planets that are lost to the 
inner disc in population synthesis models. 

7. Summary 


In this paper we have com bined previou s results and the- 


ories of planet migration (WardI 1997t iLin & Panaloizoul 

19861 Paardekooner & 

Mellema 2006h with pebble 

accre- 

tion (lormel & Klahd 

20101: iMorbidelli & Nesvornv 

20121 ; 

Lambrechts & Johansen 

120121 ILambrechts & JohansenI 

20i4 

Lambrechts et alJ 2014 

) and gas accretion (iMachida et alJ 
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Fig. 10. Categories of planets as a function of rg and tg. The different 
colours indicate the different types of planets, which are separated by 
planetary mass and their core-to-envelope mass ratio. The hot, warm, 
and cold gas giants are separated by their semi-major axis (see Table[T]l. 
The black region contains small objects with M < O.IMe. 



achieve a certain planetary mass and orbital distance to the star, 
which we used to define some distinct planetary types. 

The initial planetary seeds grow v ia pebble accretion unti l 
they reach their pebble isolation mass (iLambrechts et al.ll2014l) . 
which is when gas accretion can start. This initial growth stage 
depends on the metallicity of pebbles Z in the disc, i.e. the num¬ 
ber of pebbles available to be accreted, and on the location in 
the disc, because the pebble isolation mass depends on the as¬ 
pect ratio of the disc (eg. ITbll. T he latter changes in time as the 
disc evolves (iBitsch et al.ll2015h . As the aspect ratio of the disc 
decreases in time, the final mass of the cores become lower for 
planets that start their evolution in the disc at a later time tQ. 
Likewise, the mass of the core increases for planets that start to 
form at greater orbital distances ro, because of the flaring struc¬ 
ture of the disc, where H/r increases with increasing r and thus 
increases the pebble isolation mass. 

After the planet reaches pebble isolation mass, it can start 
to accrete gas. But before runaway gas accretion can start, the 
gaseous envelope has to contract to reach < Menv This con¬ 
traction time depends critically on the opacity inside the plan¬ 
etary atmosphere, which is not very well constrained, therefore 
future work should aim to study gas accretion and envelope con¬ 
traction of planets with gaseous atmospheres in a much better 
way. After the envelope contracted, runaway gas accretion can 
start and the planets can become gas giants. 

During the growth process of the planets, they migrate 
through the disc. While they accrete pebbles, the planets are 
small and do not significantly disturb the surroundings of the 
disc, meaning they migrate in type-I migration. This migration 
is inwards in the outer regions of the disc and can be outwards in 
the inner regions of the disc. However, these regions of outward 
migration change as the disc evolves and accretes onto the star. 
When the planet starts to accrete gas, it can also outgrow the re¬ 
gions of outward migration, and the planet moves fast inwards 
towards the star. But, when the planet becomes massive enough, 
it will start to open a gap in the disc and migrate in the slow type- 
II migration, which can keep it from reaching the inner edge of 
the disc. 

We can roughly categorise the planets by their their final 
mass Mp and orbital distance rf in the disc. The formation re¬ 
gions of these classes are displayed in Fig. fTOl 
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Hot gas giants form early in the inner regions of the disc 
(small to and small ro), grow quickly to the pebble isolation 
mass, and can then be trapped in a region of outward migra¬ 
tion, where their envelope contracts until they can accrete gas in 
a runaway fashion. At that point they have outgrown the region 
of outward migration and migrate into the inner disc, where they 
get stranded close to the central star. 

Warm and cold gas giants have a much wider parameter 
space in ro - fo to form than hot gas giants. They can form early 
far out in the disc, while they can still form later in the disc (up 
to to ^ 2.2 Myr) at moderate distances (ro ~ 15 AU). They un¬ 
dergo the same kind of growth mechanisms as hot gas giants, but 
because they form initially at a greater orbital distance or a later 
time, they need to migrate more to end up in the inner disc as a 
hot gas giant. Their larger initial semi-major axis or late forma¬ 
tion time keeps them far out in the disc. 

Ice giants are planets located in the dark blue region in 
Fig. (To] They have Me > Menv and a total mass that is higher 
than 2 Earth masses. They can form at all orbital distances, de¬ 
pending on their formation time. The formation in the outer disc 
allows them to have a larger core, because of the large pebble 
isolation mass, and the late formation time hinders an efficient 
contraction of the envelope, so they do not reach runaway gas 
accretion. In the inner parts of the disc, the pebble isolation mass 
is lower (low ///r), which prolongs the contraction of the enve¬ 
lope, resulting in a lot of planets with > Menv, which are ice 
giants by our definition. 

Ice planets are small planets that have formed in the cold 
parts of the disc rf > riee. These planets are very common in the 
outer regions of the disc, if the growth of planets starts in the 
late stages of the disc evolution (to > 2 Myr). These planets did 
not reach pebble isolation mass and have therefore only a very 
minimal gaseous envelope, which should allow them to have a 
density similar to ice giants. 

In contrast to previous simulations of the formation of giant 
planets by planetesimal accretion, we do not make any simplifi¬ 
cations regarding planet migration, nor do we have to assume an 
unreasonably high amount of solid s in the protoplanetary disc to 
form planets (iPollack et al.lll99^ . Additionally, our protoplan¬ 
etary disc evolves in time and changes not only its total mass, 
but also its temperature and density profiles as the accretion rate 
decreases. Therefore, the pebble accretion scenario suggests that 
different types of planets emerge as a result of different forma¬ 
tion time and location in the disc. 

We have shown here that pebble accretion can overcome 
many of the challenges in the formation of ice and gas giants 
in evolving protoplanetary discs. Gas giants in our model do not 
form in situ, but migrate over several AU during their formation 
process, which requires a formation of the planetary seeds far 
out in the disc. In contrast to that, ice planets and ice giants are 
more likely to form in situ compared to giant planets. In fact, 
ice planets and ice giants only form in the late stages of the disc 
evolution or in the outer parts of the disc. A late formation time 
thus gives rise to a wide variety of planetary types, akin to those 
found in our solar system, as well as in extrasolar planetary sys¬ 
tems. 
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Appendix A: Dispersal time of the disc 

Observations have shown th at the lifetime of pro toplanetary 
discs is a around a few Myr (iHartmann et alJlT998h . However, 
the fracti on of stars with discs reduces greatly with increas¬ 
ing time (lMamaiekll2009l) . In our nominal model, the lifetime 
of the disc is fixed to 3 Myr, where we follow the disc evolu¬ 
tion of IHartmann et al.l (Il998h . At 3 Myr the disc has reached a 
stellar accretion rate of M = 2 x 10“^Mo/yr. We then assume 
that the disc gets cleared by photoevaporation, so we stop our 
simulations at that time. However, the effects of photoevapora¬ 
tion a re not clear and not constrained perfectly dAlexander et al.l 
HoH, so that we now make the assumption that the disc will 
live longer, up to 5 Myr, whe re the disc will reach a n accretion 
rate of M = lx 10“®Mo/yr (IHartmann et al.l [19^ . The addi¬ 
tional lifetime of 2 Myr therefore only models the decay of the 
accretion rate from M -2x lO^^Mo/yr to M - lx 10“^Mo/yr. 

Starting from the same initial configuration as for Fig. |2] 
namely starting planets at different orbital distances when the 
disc has already evolved to fo =2 Myr, results in Fig. IA.il where 
we present the growth tracks of planets in discs that live 3 Myr 
and 5 Myr. The planets are inserted at ?£> = 2 Myr, so they un¬ 
dergo 1 Myr or 3 Myr of evolution, respectively. The evolution 
of the planets in the disc that lives longer is identical to the disc 
that lives for a shorter time, because both discs have the same 
properties. 

In the disc that lives longer, planets have more time for their 
evolution. This means that the planets, which have not con¬ 
tracted their gaseous envelope completely, can now do so and 
thus start rapid gas accretion. This is clearly seen for the planet 
in the inner system (Fig. lA.lb . which classified as an ice gi¬ 
ant in the 3 Myr disc, while it is a hot gas giant in the disc 
that lives for 5 Myr. Similar fates are shared by the planets 
that form in the outer system. Those who already started rapid 
gas accretion accrete so much gas that the planet becomes mas¬ 
sive enough that type-II migration shifts from disc-dominated to 
planet-dominated (eq.|22]), reducing its migration speed. The ice 
giants in the outer disc now have more time to accrete pebbles 
and reach their isolation mass, so that they can contract their en¬ 
velope and undergo runaway gas accretion, transforming them 
into gas giants. Even the planetary seed that would only grow to 
become an ice planet in a disc that lives 3 Myr (Fig. [2]i evolves 
into a gas giant. This also indicates that the lifetime of the disc 
does not only play any role in terminating the growth of gas gi¬ 
ants, but for planets in general. 

In Fig. IA.2I we present the ro-fo map for a disc that under¬ 
goes 5 Myr of evolution. The general shape of the different re¬ 
gions where different planetary types emerge (e.g. gas giants, ice 
giants) is very similar to Fig. |4] where the disc lived for 3 Myr 
in total. The difference is now that the formation of all plane¬ 
tary types can start at 3.5 Myr, compared to 2 Myr in the disc 
that lives 3 Myr. Otherwise, the general results do not change. 
A longer disc lifetime transforms smaller planets to larger ons, 
because they have more time to grow, and reduces the final semi¬ 
major axis of those planets, because they have more time to mi¬ 
grate. 

Appendix B: Constraints on the opacity during 
envelope contraction 

After the core reaches pebble isolation mass, it hosts a gaseous 
envelope that first undergoes contraction on a long time scale, 
before runaway gas accretion can start. During this contraction 
phase, the envelope grows slowly until M^ = Menv, which is 
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Fig. A.l. Evolution of planets that start at to = 2 Myr, but in discs that 
live up to to = 3 Myr or to = 5 Myr. The black circular symbols mark 
to = 3 Myr, which corresponds to the final positions in Fig.|^ The black 
squares mark the final position at ?o = 5 Myr. The evolution between 
3 Myr and 5 Myr is indicated by the dashed lines. Clearly planets in 
the disc with longer lifetimes continue their evolution, so that the final 
dissipation time of the disc is very important to set the final mass and 
orbital position of planets. 
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Fig. B.l. Mass of the planet as a function of time. The mass of the core 
(red line) increases via pebble accretion until isolation mass is reached 
at 350 kyr. Then the envelope starts to contract (dashed lines), which in 
increases the total mass (solid lines). The contraction time of the enve¬ 
lope depends on the opacity of the envelope where lower opacity 
results in a shorter contraction time, so that when rapid 

gas accretion can start. The planet is placed in a disc at fo = 2 Myr at 
ro = 10 AU. 
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Fig. A.2. Final masses of planets as a function of initial radius ro and 
initial time to in the disc, where Z = 1.0%, and a total disc lifetime of 5 
Myr. The lines inside the plots correspond to the same meaning as the 
lines in Fig.|4l Clearly the formation of small planets is delayed to times 
after to = 3.5 Myr. 



when rapid gas accretion can start. The accretion rate during the 
contraction phase is given by eq. [17] The opacity /Cenv inside the 
envelope is crucial in determining the accretion rate during en¬ 
velope contraction. Unfortunately, the opacity inside the plan¬ 
ets envelope is poorly constrained (llkoma et alJl200Cll) . Not only 
is the size of the grains important for determining the opacity, 
but also the temperature in the planetary atmospher e. As grain 
growt h can be efficient in planetary atmospheres (iMordasinil 
1201 4l) and the temperature of a planet accreting gas is higher than 
the surrounding disc, the opacity of the planetary atmosphere is 
reduced compared to the opacity in the disc. Additionally, the 
opacity depends on the underlying chemical composition of the 
grains, which depend on the composition of the disc and the tem¬ 
perature in the disc. 

The nominal opacity in eq. [T7|is set to A-env = 0.05 cm^/g. 
Further reductions in the opacity of the envelope result in faster 
contraction times of the envelope, so that Me < Menv is reached 


at shorter times, allowing runaway gas accretion during the 
disc’s lifetime. In fact, a close-to- zero opacity allows ra pid gas 
accretion onto cores of just IMe (iHori & Ikomall2()Tol) . When 
setting /fenv = 0.05 cm^/g, the contraction time of the envelope is 
roughly of the as the building time of the core via the pebble ac¬ 
cretion shown in Fig. ED Th is value is very close to the esti mate 
of opacity in the envelope of iMovshovitz & PodolakI (l2008h . 

In Fig. IB.ll we display the growth of a planet as a function 
of time. After the core has reached pebble isolation mass, the 
envelope starts to contract, where the opacity in the envelope /Cenv 
sets the contraction time. It clearly shows that higher opacity in 
the envelope prolongs the contraction time of the envelope. 

A prolongation of the contraction time has important con¬ 
sequences for the formation of planets, because the planet will 
then spend a longer time in type-I migration before it is massive 
enough to open a gap in the disc. This means that the distance the 
planet migrates from its initial position ro is greater for a higher 
values Menv, unless the planet is caught in a region of outward 
migration. However, this can only happen in the early evolution 
stages of the disc. In the late stages of the disc evolution, a longer 
contraction time will lead to an enlargement of parameter space 
that allows for the formation of ice giants at a few AU, because 
these planets have very small cores and therefore do not contract 
an envelope. 


Appendix C: Planet formation in the MMSN 

To emphasize the importance of a realistic disc structure for the 
formation of planets, we tested our planet formation model also 
with a MMSN disc (IWeidenschillind [1 9771: Havashil 19 811). A 
simple power law disc was also used in Lambrechts & .lohansenl 
(l2014h . but with slightly different parameters compared to the 
MMSN. Here we directly probe the MMSN, where the discs sur¬ 
face density is given by 





(C.l) 
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Fig. C.l. Final masses of planets as a function of initial radius ro and 
initial time to in the MMSN with a lifetime of 3 Myr. The lines inside 
the plots correspond to the same meaning as the lines in Fig.|4] We again 
us e a metallicity of pe bbles of Z = 1.0%. In contrast to the disc model 
of iBitsch et all i20lk) . the MMSN disc allows 0.1 AU < rf < 5.0 AU 
only for a very small amount of parameter space. Either the planets fall 
close to the central star (rf <0.1 AU) or they stay outside of 5 AU. On 
the other hand, the parameter region that allows for the formation of ice 
giants is greatly increased compared to Fig.|4] 

The aspect ratio of the disc follows 
H ( r 

- = 0.033 (-) , (C.2, 

which indicates that there are no planet traps due to the entropy 
driven corotation torque, b ecause those traps exist only when 
H/r decreases with radius (iBitsch et al.ll2014l) . The time evolu¬ 
tion of the surface density is given by 

^2 = 2oexp(—] , (C.3) 

\ ^disc / 

where raise = 3 Myr and Sq = 1700g/cm^. From this disc struc¬ 
ture we model the formation of planets, and the resulting ro-fo 
diagram is shown in Fig. ED 

The resulting planets look dramatically different from those 
in Fig. m In the inner regions of the disc (rp < 10 AU), the 
planets do not reach rapid gas accretion, because the migration 
in the inner disc is so rapid (because of the high surface density) 
that the planets arrive at rf < 0.1 AU before Me < Menv Planets 
that form farther outside can reach this stage (Me < Menv), but 
their migration is still so rapid that they end up at rf < 0.1 AU. 
Nearly a third of the ro-fo parameter space results in planets with 
rf < 0.1 AU. Only a very small band allows for the formation 
of planets with 0.1 AU < rf < 5.0 AU, where most of these 
planets have Me > Menv, indicating that the formation of ice 
giants outside of 1 AU and inside of 5 AU is very hard in this 
disc model. 

However, outside of 5 AU, a very large band allows for the 
formation of ice giants (green region in Fig. 1C. II) . These ice gi¬ 
ants have a much larger number of solids (~ 20 - 30Me) com¬ 
pared to Fig. m This is caused by the higher aspec t ratio in the 
MMS N disc, compared to our nominal disc model (IBitsch et al.l 
1201 5l) . which increases the pebble isolation mass (eq. [Tbl l. Ad¬ 
ditionally, in the MMSN disc model, only the surface density 
evolves in time, while the aspect ratio does not evolve, allow¬ 
ing a high pebble isolation mass in all disc evolution stages. In 
combination with the efficient pebble accretion, planets in this 
disc model will have a higher core mass than in the nominal disc 
model used in Fig.|4] 


Appendix D: The heating torque 

A recent study by iBenftez-Llambav et alJ (12015l) shows that a 
small planet (Mp < 5Me) that accretes very quickly (mass dou¬ 
bling times shorter than 60 kyr) alters the structure of the sur¬ 
rounding gas disc in such a way that it can migrate outwards 
in the disc. The accreting material onto the planet causes asym¬ 
metries in the temperature of the disc close to the planet that 
produce a force that counteracts inward migration. This effect 
concerns planets that are so small that they had not yet reached 
the region of outward migration (see Fig. [T]). Even if the mass 
doubling time is large, e.g. 300 kyr, the heating torque has an ef¬ 
fect, because it reduces the inward speed of the planet compared 
to the nominal type-I migration rate. 

This so-called heating torque is a strong function of the ac¬ 
cretion rate onto the planet (only mass doubling times shorter 
than 60 kyr lead to outward migration) and of the opacity 
of the disc, which determines the cooling in the surround¬ 
ing disc. Higher opacity suppresses cooling, which leads to a 
stronger effect of the heating torque. In their standard set-up, 
iBenftez-LIambav et al.l (l20I5l) use a constant opacity of I cm^ /g, 
and the planet has a mass of 3 Me. The effect of the heating 
torque also depends on the mass of the planet. Unfortunately, 
this new effect is not quantified in a function that contains all rel¬ 
evant parameters (planetary mass, opacity of the disc, and mass 
doubling time). We therefore make a simple test, where we as¬ 
sume that the heating torque is just a function of the accretion 
rate and do not investigate effects of opacity and planetary mass 
in a detailed way. For that we use the data in the methods section 
of iBenItez-Llambav et al.l (l2015l) . 

We assume that the effect of the heating torque is equally 
strong until the planet reaches 5 Me, after which the heating 
torque is not relevant any more. Additionally, we assume just for 
the calculation of the heating torque that the opacity is equal to 1 
cm^/g, which is actually not the case in our disc. In fact, only in 
the inner parts of the disc, so close to the ice line, is the opacity 
higher than 1 cm^/g in our disc model. Depending on the initial 
semi-major axis ro and the initial time fo, the mass doubling time 
in our simulations can be a few kyr (small planet with ~ 0.5 Me, 
small ro, and small fo) up to a few 100 kyr (large planet, larger 
ro, and larger fo). 

In Fig. ID. II we display the evolution tracks of planets with 
and without the heating torque. In the example where the planets 
start at 5 AU, they are also inserted into the disc at an early time 
(fo = 100 kyr). In the early disc, the region of outward migra¬ 
tion is still very large (Fig. [TJ, so that planets that are massive 
enough (a few Earth masses) to reach this region and can grow 
there without migrating inwards. In this example, because ro and 
fo are small, the heating torque easily generates outward migra¬ 
tion, and the planet ends up in the region of outward migration. 
However, the planet also reaches the region of outward migra¬ 
tion when the heating torque is not taken into account, because 
the growth by pebble accretion is so fast that the planet does not 
migrate much during this time. Both planets then end up on the 
same evolution track, because they reach the region of outward 
migration where their previous migration history no longer mat¬ 
ter. 

The planets starting at 20 AU in Eig. IDTI also start later in 
the disc at fo = 2 Myr. There the planets do not reach a region 
of outward migration. Additionally, the mass doubling time be¬ 
comes longer than 60 kyr for most of the growth stages of the 
planets, so that the heating torque does not prevent inward mi¬ 
gration, but just slows it down. However, even this slowing down 
has an effect on the final orbital mass and position of the planet. 
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Fig. D.l. Two examples of the evolution of planets with and without 
heating torque. Both planets starting at 5 AU start at an early initial time 
to = 100 kyr, while the planets starting at 15 AU start at a late initial time 
of fo = 2 Myr. The small black dots indicate a time difference of 200, 
400, 600, and 800 kyr from the starting time. The big black dots indicate 
an evolution time of 1 Myr. The planets starting early and in the inner 
disc only evolve for ~ 500 kyr total. With the fast growth rates of pebble 
accretion, it seems that the heating torque only has minimal effects on 
the final configuration of evolving planets. 


The planet that grows when we take the heating torque into ac¬ 
count has a somewhat larger final semi major axis, but reaches a 
similar final planetary mass. The reason for this effect not seem¬ 
ing too important here is that small mass bodies migrate quite 
slowly, and they also grow quite quickly thanks to pebble accre¬ 
tion, so that they reach 5 Me rather quickly, so that the planet 
experiences the effects of the heating torque for only a small 
amount of time. 

The general results of our simulations, for example Fig. H] 
still holds, even when taking the heating torque into account. The 
only thing that changes is that the final orbital positions rf of the 
planets are slightly larger, meaning that a smaller ro is needed to 
achieve a certain final orbital position rf. In the simple test here, 
the initial semi-major axis can be up to 2 AU smaller in order 
to achieve the same rj for a given fo compared to simulations 
where the heating torque is not taken into account, making the 
effect in a global picture not that important here. 

This effect seems to not be too important in the pebble accre¬ 
tion scenario because of the fast growth rates. The fast growth 
rates for reaching Mp > 5 Me mean that there is only a short 
time for the planet to actually spend migrating. The resulting 
difference between simulations with and without heating torque 
where the accretion rates are high is therefore minimal. How¬ 
ever, the heating torque could have substantial effects when the 
accretion rate is low, and the heating torque redu ces the torque 
causing inward migration. In ICossou et al.l (l2014t) the growth of 
planetary cores is modelled in N-body simulations via the colli¬ 
sion of planetesimals. There the production of the cores of giant 
planets at large orbits is hindered by the fact that planets mi¬ 
grate inwards too fast below the region of outward migration 
(because the planets are not massive enough and do not grow 
fast enough), which results in many super-Earth types of planets 
instead of giant planets. Including the heating torque this picture 
could change, because the inward migration speed is reduced for 
these small objects, which could allow them to grow more and 
reach the region of outward migration. Additionally, if the metal- 
licity caused by pm sized dust grains is greater than the 0.5% as 
in our work, the heating torque could become more prominent, 
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Fig. E.l. Final masses of planets as a function of initial radius ro and 
initial time fo in the disc with Z = 1.0% (top) and Z = 1.5% (bottom). 
The lines inside the plots correspond to the same meaning as the lines 
in Fig. |4] Plan ets below the purple line are subject to dynamical coro¬ 
tation torques jPaardekoopenUof^. while planets below the g reen line 
are subject to type-III mieration jMasset & Papaloizoull200^ . indicat¬ 
ing that their final orbital position does not correspond exactly to the 
position indicated by the black lines. 


because it strongly increases with increasing opacity, which di¬ 
rectly scales with the metallicity. 


Appendix E: Type-Ill migration and dynamical 
corotation torque 

When a planet is big enough to carve a gap around its orbit, the 
horseshoe region becomes depleted. However, when the mass 
pushed away by the planet is higher than the mass of the planet 
itself (the co-orbital mass deficit), the planet can undergo a rapid 
change of semi-major axis, which is described as type-III mi¬ 
gration (iMasset & Papaloizoull2003l) . Migrating planets also ex¬ 
perience dynamical torques, which are proportional to the mi- 
gration rate and dep end on the background vortensity gradient 
(|PaardekooDerl20 fil) . Here we want to discuss how these effects 
influence our models. The description of how these mechanisms 
work can be found in section l24l 

In Fig. IE. 1 1 we have indicated the planets that would either 
undergo type-III migration or are subject to dynamical coro¬ 
tation torques (everything below the green and purple lines in 
Fig. IE.lb resulting in the fact that the final orbital positions in 
these regions are not exactly marked by the black lines. The final 
orbital positions are most likely closer to the host star, because 
the runaway type-III migration follows the direction of motion, 
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which is inwards in regions outside of r > 12 AU (Fig.[B. This 
process only occurs in the early stages of the disc evolution, 
where the disc is still massive. In the later stages of the disc 
evolution, the disc is less massive, and these effects are not sig¬ 
nificant any more, making the predictions of our model reliable 
in this part of parameter space. 

Simulations with a larger number of pebbles can form gas gi¬ 
ants more easily and in a larger parameter space in ro-fo (Fig-lS. 
A larger number of pebbles results in a faster growth of the plan¬ 
etary core, which then can contract its envelope faster as well, 
leading to a total faster growth of the planet. This means the 
planet will reach the mass where it starts to open a gap partially 
in the disc earlier as well. However, in an earlier stage of evolu¬ 
tion, the disc is also more massiv, indicating that the co-orbital 
mass deficit is larger as well, making the planet more prone to 
type-III migration. The same effect applies for the dynamical 
corotation torque. In the situation of a larger number of pebbles 
in the disc, a greater number of planets in the ro-fo parameter 
space are affected by type-III migration and dynamical corota¬ 
tion torques (Fig. ED. 
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